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Abstract. We make the first steps towards a generic theory for energy spreading 
and quantum dissipation. The Wall formula for the calculation of friction 
»v/l , in nuclear physics and the Drude formula for the calculation of conductivity 



in mesoscopic physics can be regarded as two special results of the general 
formulation. We assume a time-dependent Hamiltonian H{Q, P;x{t)) with 
x{t) = Vt, where V is slow in a classical sense. The rate-of-change V is not 
necessarily slow in the quantum-mechanical sense. The dynamical variables 
(Q, P) may represent some 'bath' which is being parametrically driven by x. 
, This 'bath' may consist of just few degrees-of-freedom, but it is assumed to be 

^ti* . classically chaotic. In case of either the Wall or Drude formula, the dynamical 

variables (Q, P) may represent a single particle. In any case, dissipation means 
Cn ' an irreversible systematic growth of the (average) energy. It is associated with 

^_^ ' the stochastic spreading of energy across levels. The latter can be characterized 

Cj> ^ by a transition probability kernel Pt{n\m) where n and m are level indices. This 

(J"^ , kernel is the main object of the present study. In the classical limit, due to 

the (assumed) chaotic nature of the dynamics, the second moment of Pt{n\m) 
j^ ' exhibits a crossover from ballistic to diffusive behavior. In order to capture this 

e' crossover within quantum- mechanics, a proper theory for the quantal Pt{n\m) 

should be constructed. We define the V regimes where either perturbation 
theory or semiclassical considerations are applicable in order to establish this 
crossover. In the limit ?i — ► perturbation theory does not apply but semiclassical 
considerations can be used in order to argue that there is detailed correspondence, 
during the crossover time, between the quantal and the classical Pt{n\m). In the 
perturbative regime there is a lack of such correspondence. Namely, Pt{n\m) 
is characterized by a perturbative core-tail structure that persists during the 
crossover time. In spite of this lack of (detailed) correspondence there may 
jr^ 1 be still a restricted correspondence as far as the second-moment is concerned. 

5^ ' Such restricted correspondence is essential in order to establish the universal 

C^ ' fluctuation-dissipation relation. 



1. Introduction 

1.1. Definition of the problem 

We consider in this paper a system that is described by an Hamiltonian H{Q, P; x) 
where (Q, P) are canonical variables and a: is a parameter. It is assumed that 
T-i{Q,P;x) with x = const generates classically chaotic motion. We are mainly 
interested in the case of time dependent x(t). However, it is assumed that i = 1/ is a 
classically small velocity. The notion of classical slowness will be defined in SecJS. The 
theory that we are going to present is quite general. In some particular applications 
x{t) may represent, for example, a time-dependent electric field. However, the theory 
is best illustrated by considering the 'piston' example: In this example x represent the 
position of a small rigid body that is translated inside a large cavity, and (Q, P) are 
the coordinates of a tiny gas particle. See Fig.|l|. 
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Figure 1. The 'piston' example: Tlie slow degree of freedom is the 'piston', and 
the 'bath' consists of one gas particle. In this paper the position of the 'piston' 
x(t) = Vt is treated as a classical parameter. Dissipation means a systematic 
growth of the bath-energy. For simplicity we assume throughout the paper that 
the conservative work is zero. This is not the case in the right illustration. 



It is assumed that initially the system is characterized by some energy distribution 
p{E). In particular we may assume a microcanonical preparation. For V = energy 
is a constant of the motion, and therefore the energy distribution p{E) will not change 
as a function of time. On the other hand, for V^ 7^ the energy will be re-distributed 
and p{E) will become time dependent. In this paper we are interested in the study of 
this time dependence. Of particular interest is the time dependence of the first and of 
the second moments. A systematic increase of the average energy has, by definition, 
the meaning of dissipation. In case of the 'piston' example, dissipation means that 
the gas particle is being 'heated up'. 



1.2. Restrictive sense of 'Quantum dissipation' 

The subject of this paper is the quantum-mechanical (QM) theory of energy-spreading 
and dissipation, as defined in the previous subsection. In short we may say that we are 
interested in the theory of Quantum Dissipation. However, it is important to realize 
that we are using the term 'Quantum Dissipation' in a quite restrictive sense. This is 
because of mainly two reasons: (a) We assume a classical driving force; (b) We are not 
considering a many-body bath. Note that an infinite number of degrees-of-freedom is 
not important for having stochastic behavior: this is the main idea behind the term 
'chaos' when applied to dynamical systems. We can have dissipation even if {Q,P) 
represent a few degrees-of-freedom 'bath'. 

The interest in Quantum Dissipation is very old pi 0, S, 0, 0, pi Q . However, in 
most of the literature, the term 'Quantum Dissipation' is used in a more general sense. 
Namely, x becomes a dynamical variable, and one looks for its reduced dynamics. 
Thus, in most of the literature, dissipation-of-energy becomes only one aspect of a 
much more complicated problem. The 'grand problem' includes, besides 'dissipation', 
other issues such as 'dephasing' and 'thermalization'. It also should be noticed that 
the standard literature usually adopts an effective-bath approach (see subsection 1.4) 
or other effective formulations [g| that do not necessarily reflect the actual dynamics 
of the bath degrees-of-freedoms. Important exceptions are works such as Q and |g^ . 

Of particular interest is the 'piston' model. If the 'piston' is treated as a dynamical 
object, than its reduced dynamics is called 'quantal Brownian motion' (QBM). 
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According to our (restricted) definition, 'dissipation' means systematic irreversible 
growth of the bath-energy. In case of an un-driven Brownian particle, the 'dissipation' 
is balanced eventually by 'noise' leading to 'thermalization'. In the QM case the 
issue of 'irreversibility' is more complicated because we may have 'recurrences'. The 
relevant time scale for these recurrences is the Heisenberg time for the combined 
BrownianParticle-GasParticle system. This latter time scale may be extremely large 
if the Brownian particle has a large mass. 

In this paper x is not a dynamical degree of freedom, and therefore the 
'recurrences' that have been mentioned in the previous paragraph are not an issue. 
(It is as if we assume that the 'piston' has an infinite mass, hence the frequent use 
of the term 'moving walls'). For V — the Hamiltonian is time- independent, and 
we will have recurrences that are associated with the dynamics of the GasParticle 
(alone). The remnant of this latter type of recurrences is QM-adiabaticity, which we 
are going to discuss soon. Another type of 'recurrences' are associated with periodic 
driving and are discussed in Sec.^. It should become clear from the above discussion, 
and subsection 1.6 below, that 'recurrences' are not an important issue in this paper. 

1.3. The classical theory of dissipation 

The classical understanding of the dissipation-process is based mainly on the works 
of HI 0' H H' Q ^^"^ followers. We are going to sketch briefly the main idea of 
the classical theory, and the associated derivation of the fluctuation-dissipation {J-T>) 
relation. In the time-independent case {V = 0) the motion of {Q(t), P(t)) is irregular 
due to the chaotic nature of the dynamics. We shall denote the ergodic time by t^^g. 
We can define a fluctuating quantity J-'{t) = —{dTi/dx) that has stochastic features. 
The intensity of these fluctuations will be denoted by v. In the classical case T{t) is 
essentially like noise whose correlation time r^i is smaller than or equal to ie^g. In the 
time- dependent case [V ^ 0) energy is not a constant of the motion and consequently 
the energy distribution p{E) becomes time dependent. It is argued that for t ^ i„g 
the energy distribution satisfy a diffusion equation. The energy-dependent diffusion 
coefficient will be denoted by De- It turns out that quite generally D^ = \vV^ . 
Associated with this diffusion is a systematic growth of the average energy. This 
systematic growth of energy is due to the £'-dependence of the diffusion process. The 
rate of energy growth will be denoted by Q. It can be written as Q = jiV^ . 

The considerations above lead to the conclusion that in the classical case the 
dissipation is of ohmic nature (Q ex V"^). The dissipation coefficient is denoted by 
/i. It is implied that the fluctuating quantity !F{t) has a non-zero average, namely 
(JF) = —fiV. In the 'piston' example the latter represents the 'friction' force that 
is experienced by the moving object. The considerations above also imply that the 
analysis of dissipation is reduced to the study of energy spreading. The difficult issue 
is to establish a stochastic energy spreading with a coefficient D^ — \vV^ . Then, the 
!FD relation between /i and the noise intensity v follows as an immediate consequence. 
If p{E) is a canonical distribution (which is not necessarily the case) then the !FD 
relation reduces to the familiar form p, = v / {2k^T) where T is the temperature. 

1.4. The effective-bath approach to Quantum Dissipation 

The most popular approach to 'Quantum Dissipation' is the effective-bath approach 
[El, 0, pi, 1^. When applied to 'our' problem (as defined in the first subsection) it means 
that the chaotic (Q, P) degrees-of- freedom are replace by an effective-bath that has 
the same spectral-properties. This may be either harmonic-bath (with infinitely many 
oscillators) or random-matrix-theory (RMT) bath [pO[ . 
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Generic classical parameters (t„, , ly) 
T^i = classical correlation time. 
1^ = intensity of fluctuations 

Generic quantal parameters (A, b, a, h) 

A = mean level spacing of the eigen-energies {En} 

b = Dimensionless bandwidth of the matrix (dTL/ dx)nm 

a = Root-mean-square of in-band matrix elements of [dT-L/dx) 

Semiclassical relations 

T,, = 27rn/(6A) 

V 



(27r?i/A) cr2 



Linear response theory 

De = \yV^ = {TTh/A) a^V^ 

Ohmic dissipation 

d{n)/dt = ^y2 
/i = jp[i^] 

Primary dimensionless parameters {b, Upr) 



= (l/h) 0^ V = b-"^ (27r?i/A)2 {a/h) V 



Table 1. Overview of the common theory for dissipation. Two generic parameters 
should be specified for the classical theory, while four are required for the QM 
theory. Note that V = x always appears in the combination uV^ or aV , and 
therefore it should not be counted as an additional (independent) parameter. 
The two classical parameters can be expressed in terms of the QM parameters via 
semiclassical relations. In the absence of well defined classical limit (as in the case 
of RMT models) this relations can be regarded as definitions. The so called Kubo- 
Greenwood result of linear response theory can be obtained using FGR picture, 
and it coincides with the classical expression. General considerations lead to a 
fluctuation-dissipation (.?2?) relation between /i and v. An important observation 
of this paper is that the validity of the linear-response approach is controlled by 
the dimensionless parameter Vp^.- 



It turns out that quantal-classical correspondence (QCC) is a natural consequence 
of this procedure: The dissipation coefhcient ji turns out to be the same classically and 
quantum-mechanically. In order to explain this point let us use the Caldeira-Leggett 
notations 0. The distribution of the frequencies of the bath-oscillators is characterized 
by an ohmic spectral-function J{'jj) — rjuj. The classical analysis leads to a friction 
force with a coefficient /i = 77, and white noise whose intensity is j^ = 2rik^T . Using 
Feynman- Vernon ||l| formalism one obtains the same value jjl — rj in the QM case. 
The quantal noise is characterized by an ?i-dependent power-spectrum, but the noise 
intensity v is defined as the w = component, and it is still equal to 2rjk^T . Hence 
the classical Jl^ relation /i = v/{2kjiT) holds also in the QM case. 

The effective-bath approach will not be adopted in this letter since its applicability 
is a matter of conjecture. In this paper we want to have a direct understanding of 
quantum-dissipation. 



Theory of Quantum Dissipation 5 

1.5. The QM theory of Dissipation 

Quantum- mechanics introduces additional energy scales, as well as additional 
parametric scales into the problem (See Table 2). Consequently there are few V 
regimes in the QM theory (See Table 3). The QM-adiabatic regime is quite well 
understood. We shall discuss this regime only briefly since it is not related to the 
main concern of this paper. The further distinction between the QM-slow regime and 
the QM-fast regime is the main issue of this paper. 

Let us assume that initially the energy is concentrated in one particular level. For 
extremely slow velocities (ulz ^ 1) and relatively long time the energy will remain 
mainly concentrated in the initial level. This is the QM-adiabatic approximation. The 
term 'QM-adiabaticity' is a beat confusing, because it actually does not correspond 
(in the h ^ limit) to adiabaticity in the classical sense. Maybe a better term 
would be 'perturbativc localization'. In the QM-adiabatic regime Landau-Zener 
transitions between neighboring levels constitute the predominant mechanism for 
energy spreading. This mechanism does not correspond to the classical mechanism of 
energy-spreading. The QM-adiabatic regime is a genuine quantal regime. 

For higher velocities (wlz ^ 1) it is essential to take into account transitions 
between non-neighboring levels. The contribution of near- neighbor transitions to the 
energy spreading becomes negligible rather than predominant. An obvious approach 
for the study of energy spreading would be to adopt a Fermi-golden-rule (FGR) picture. 
FGR is one possible picture of perturbation theory. The same results for Z?e and fj, 
can be derived by using other, equivalent formulations of perturbation theory. The 
most popular variation is known as 'linear response theory' or as 'Kubo-Greenwood 
formalism'. Whatever version of perturbation theory is being used the standard result 
is always the same (See Table 1). It should be realized that the standard result is in 
complete correspondence with the classical result, and it becomes identical with the 
classical result upon taking the formal limit ?i — > 0. 

1.6. Specific motivation for the preset study 

Reading some of the early literature one gets the impression that quantum dissipation 
is conceptually well-understood. Specifically, it looks as if the perturbativc methods 
are effective for the purpose of constructing a general theory. However, this is a wrong 
impression. A general theory of energy spreading is still lacking, a-fortiori there is 
no general theory of quantum dissipation. This point becomes most evident once we 
read the work by Wilkinson and Austin (WfcA) |10| . Their observations constitute 
the original motivation for the present study |13| ]. 

W&A ||l^ have defined two important dimensionless parameters that are 
associated with the velocity V. These are, (using our notations), the scaled velocity 
i>Lz, and the scaled velocity Vrmt- The QM-adiabatic regime is distinguished by the 
condition Wlz ^ 1, where we have the relatively simple picture of spreading due to 
Landau-Zener transitions. At higher velocities (wlz ^ 1) the QM-adiabatic nature 
of the dynamics is lost, and the Landau-Zener picture no longer apply. In order to 
extend the perturbative treatment to such higher velocities W&A have suggested to 
adopt an innocent-looking RMT assumption. As long as the velocity is sufficiently 
slow (wrmt ^ 1, but still Wlz ^ 1) a classical-like result for D^ is obtained. On the 
other hand, once v„mt ^ 1, the classical- like expression for D^ no longer holds. It is 
modified in such a way that correspondence with the classical result is being lost! 

Obviously, W&A have realized that the above conclusion is inconceivable. We 
will have to understand what is wrong with their innocent-looking RMT assumption. 
We shall argue that Wrmt '^ 1 does not mark a crossover to a non-classical regime. 
Rather, we shall find out that there is a different dimensionless parameter (wpr) that 



Theory of Quantum Dissipation 



Energy Scales: 

A (X h — mean level spacing of the eigen-energies {-En}- 

Af, — 6A = 2Trh/T^i — bandwidth of the matrix [d'H/dx)nm 
Asc oc h^'^ — semiclassical width of Wigner function. 

Parametric scales: 

Sx'j} = parametric correlation scale of the x-dependent Hamiltonian 
(Ja;^'" = (^/o') °c hP''^ •* " = The 5x required to mix neighboring levels. 
Sx^^t — v6(A/cr) oc h — The Sx to mix all the levels within the bandwidth. 
Sxsc fx h'^'^ — The dx required to get detailed QCC. 

Temporal scales: 

r^i = Classical correlation time of JF(t). 

iorg = Ergodic time of the classical chaotic motion. 

if,„ = i'/{fiV)'^ — Breaktime of the classical adiabatic approximation. 

r^" = Sx^^/V = The time it takes to mix neighboring levels . 

iprt = Ultimate breaktime of the QM perturbation theory. 

igj„ — Breaktime of the QM sudden approximation. 

tn — 27r?i/A = Time needed to resolve individual levels (Heisenberg time). 



Table 2. Various scales in the theory of energy spreading. The generic h 
dependence is indicated in most cases. The parametric scales and the temporal 
scales are associated with the kernels P{n\m) and Pt{n\m,) respectively. The 
determination of tp^t and iadn is an important issue of this paper. Their 
dependence on V is illustrated in Fig.bl It should be realized that r^i can be 
defined, from a purely QM point of view, as the time which is required in order to 
resolve the energy scale A;,. Similarly iadn is defined as the time which is required 
in order to resolve the spreading profile. The time t^^^ can be either equal or 
shorter than Td. 



controls the route towards quantal-classical correspondence (QCC). 

1.7. Main claims of this paper 

The purpose of this paper is to describe the time-evolution of the energy spreading 
in the various velocity regimes (See Table 3). The various 'scenarios' are graphically 
illustrated in Fig.H. The main claims of this paper are implied by this illustration. 
Disregarding the QM-adiabatic regime we are motivated by the following two 
questions: 

• What is the regime where FGR/RMT picture is valid? 

• What is the regime where QCC considerations are valid? 

In particular we would like to know whether there is a 'clash' between FGR/RMT 
considerations on the one hand, and QCC considerations on the other hand. 

The main object of this paper is the transition probability kernel Pt{n\m). The 
variable m denotes the initial energy-state of the system, and n stands for one of the 
instantaneous energy-states at a later time t. This kernel is well defined quantum- 
mechanically as well as classically. The energy distribution pt{E) can be obtained by 
operating with the kernel Pt{n\m) on the initial microcanonical preparation pt=o{E), 
and making a simple change of variables n i—> E. An important distinction in this 
paper is between restricted QCC and detailed QCC. Detailed QCC implies that the 
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Classical slowness condit 


ions: 




Ft,, < Sxl' 








Trivial condition 


Tci < if„ 








Non-trivial Condition 


Quantal regimes: 










Vtu < Sxl"' 








QM-adiabaticity (extremely slow velocities) 


Vt,, < 5x^,, 








QM-slow velocity (linear response regime) 


Vn, > fesc 








QM-fast velocity (semiclassical regime) 


Scaled velocities: 










Wsc = \/2-De t,, 


/Asc 






= V 1 (fesc/r.O 


WpH = VSi'E T,i 


/ Ab 






= V 1 {8x,Jt^,) 


71 — /)1/^ 71 


= T,i 


/ 


_qm 
' C 


^ v/{5xr/T.,) 


WlZ = fe'^''' VPB 


= ^H 


/ 


rr 


- vi{5xr/t^) 



Table 3. Definitions of the various V regimes in the theory of energy spreading 
and dissipation. The classical slowness condition is always assumed to be satisfied. 
In the QM case we distinguish between the regimes of QM-adiabaticity (extremely 
slow velocities), QM-slow velocities, and QM-fast velocities. Dimensionless 
(scaled) velocities can be defined in order to distinguish between the various 
regimes. For reasonably small h we have Vlz 2> i^rmt ^ i'pr 3> J^sc- In the 
classical limit all of them ^ 1. The condition Dlz <S 1 defines the QM-adiabatic 
regime. The condition Dpr <g 1 define the regime of QM-slow velocities. The 
condition wsc ^ 1 defines the regime of QM-fast velocities. The parameter Urmt 
has been introduced in |10[ , and we are going to explain that it determines the 
limitation of an over-simplified RMT approach. 



quantal Pt{n\m) is similar to the classical Pt{n\m). We shall see that detailed QCC can 
be established in an intermediate time regime provided the velocity V is large enough. 
In the absence of detailed QCC we still may have restricted QCC. The latter implies 
that only the first and the second moments of the corresponding distributions (quantal 
versus classical) are similar. Restricted QCC is sufficient in order to guarantee QCC 
as far as the diffusion coefficient D^ and the dissipation coefficient ^ are concerned. 
Our main statements are: 

• The FCR picture implies restricted rather than detailed QCC. 

• The FCR picture is valid in the regime Wp^ <^ f . 

• Detailed QCC considerations are valid in the regime Vsc ^ f • 

It should be realized that the FCR picture is not valid in the regime Wpr 3> 1. However, 
this does not necessarily imply that the standard Kubo-Greenwood result is not correct 
there. On the contrary: In the the limit ?i ^ we have detailed QCC {vsc ^ !)> and 
at the same time Kubo-Greenwood result simply coincides with the classical result. 
Thus we may say that for z;sc ^ 1, the standard Kubo-Greenwood result is not valid 
but correct. The distinction between 'valid' and 'correct' is crucial here: A correct 
result sometimes follows from using wrong assumptions. In the intermediate regime 
(iipR S> 1 but Wsc «C 1) neither FCR nor QCC consideration apply and we may have 
qualitatively different results for /i. It is suspected |2j, but not yet proved in the 
present context, that some artificial RMT models, that does not possess a well defined 
classical limit, may exhibit for v^^ ^ 1 a significantly different behavior compared 
with the expected FGR or classical result. This latter observation is in the spirit of 
[pO[, but it is quite different as far as details are concerned. 
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It is important to understand what is the origin of the FGR picture vahdity 
condition Wpj, <ti 1, and why it is different from the condition Wrmt ^ 1 that has 
been suggested in |10|. Again we assume that initially the energy is concentrated in 
one particular level. We shall argue that in order to determine -De it is important to 
estimate how many levels are mixed non-perturbatively at the time t ~ Td. If the 
related parametric change 6x = Vt, does not mix neighboring levels, then we are on 
"safe ground" of standard first-order perturbation theory (FOPT), and we can trust 
completely the FGR picture. Such circumstances are guaranteed by the condition 
I'rmt ^ 1. If fRMT ^ 1 we have a breakdown of the standard FOPT picture, but this 
does not imply that the FGR picture becomes non-valid. It turns out that the FGR 
result for transition rate between levels is valid on "large" energy scales, even if there 
is non-perturbative mixing of levels on "small" energy scales. This is true as long 
as the "small" scale is much smaller compared with the bandwidth Ah of first-order 
transitions. Such circumstances are guaranteed by the condition WpR <C 1. 

1.8. The 'piston' example - The wall formula 

We are using throughout this paper the 'piston' model of Fig.n^ as an illustrative 
example. It should be noticed that for simplicity of presentation we picture the 'piston' 
as a small moving obstacle whose motion is constrained to be in one space direction. 
From purely linguistic point of view 'piston' implies also hermetic closure along the 
margins. We do not assume such a closure. 

Application of the JTJ relation in order to get an expression for the dissipation 
coefficient /i leads to the 'wall' formula. This formula has been originally derived using 
kinetic considerations [^ , and only later using other approaches [^, |ll| , including the 
JT> approach that we are using here. 

In the proper classical limit (taking h-^0, while all the other parameters are 
held fixed) the walls of the 'piston' always become 'soft', meaning that De-Broglie 
wavelength becomes much smaller than the penetration distance. The hard walls 
limit (meaning that the penetration distance is taken to be zero, while h is being kept 
fixed), is non-generic. It is important to understand the consequences of taking this 
limit. In the hard wall limit (dTi./dx)nm is not a banded matrix, and the (generic) 
problem of having a non-perturbative regime for ?i — > is being avoided. The non- 
generic features of the hard wall limit are possibly responsible for some prevailing 
miss-conceptions, and in particular to the illusion that perturbative techniques can be 
used in order to get a general theory for 'quantum dissipation'. This is the reason for 
the inclusion of a quite detailed discussion of the 'piston' example. The consequences 
of taking the hard- wall limit, as well as other non- generic features of the 'piston' 
example are further discussed in the concluding section and in p4| . 

1.9. The 'mesoscopic' example - Drude formula 

Another physical examples that can be treated by the general theory of dissipation 
is taken from the realm of mesoscopic physics. Consider the case where x is the 
magnetic flux via a ring. The velocity V — x has then the meaning of electro-motive- 
force. Let us assume that the ring contains one charged particle (Q, P) that performs 
diffusive motion. Ohmic dissipation (Q = M^^) means that the charged particle gains 
kinetic energy, where the dissipation coefficient fi is just the conductivity of the ring. 
Equivalently, having {J-') = —fiV just means that the drift velocity along the ring is 
proportional to the electro-motive-force. It is a trivial exercise to get Drude formula 
from the general JT? relation. The advantage of this procedure is that the derivation 
can easily be extended to the case where the motion of the charged particle is chaotic 
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rather then diffusive. It should be noted that in actual circumstances the charged- 
particle is an electron, and its (increasing) kinetic energy is eventually transfered to 
the vibrational modes (phonons) of the ring. The latter process (that leads to Joule 
heating) is 'on top' of the generic dissipation problem that we arc going to analyze. 

1.10. Overview of the paper 

This paper divides roughly into three parts. The appendixes (A-J) should be 
considered an integral component in the reading of the main text. The reason for 
transforming some of the sections into appendixes was the desire to maintain a simple 
logical flow. We turn now to give a brief description of the paper. The first part of this 
paper (sections 2-7), in a superficial glance, looks like a review. However, inspite 
of its textbook style, it is not a review. It gives the necessary introduction for the 
later QM analysis, and in particular it contains a careful examination of the various 
assumptions that are involved in the common approaches to the theory of dissipation. 
(It turns out that a satisfactory presentation is lacking in the existing literature) . The 
main items of the first part are: 

• The crossover from ballistic to diffusive energy spreading. 

• Precise formulation of the classical slowness conditions. 

• Brief description of the derivation of the !FD relation. 

• Critical discussion of the QM linear-response theory. 

• Critical discussion of the standard FGR picture. 

• The wall formula generalized to arbitrary dimensionality (d = 2, 3...). 

The second part of this paper (sections 8-10) contains the precise formulation of the 
theory and an overview of the general picture. The main items of the second part are: 

• Definitions of the kernels P{n\m) and Pt{n\m). 

• The stochastic description of the energy spreading process. 

• Restricted QCC versus detailed QCC, and the classical approximation. 

• Overview of the dynamical scenarios in the different V regimes. 

The third part of this paper (sections 11-20) gives a detailed presentation of 
perturbation theory and RMT considerations. The main items are: 

• The Schroedinger equation in the x-dependent basis. 

• The QM-sudden approximation and parametric evolution. 

• The over-simplified RMT picture. 

• An improved version of perturbation theory. 

• The core-tail structure of the spreading kernel. 

The paper is concluded (section 21) by pointing out some important questions that 
have been left open. Some future directions for research arc indicated. 

1.11. The need for a generic theory 

Our main interest in this paper is to construct a generic theory for energy spreading 
and quantum dissipation. In particular we want to define the conditions for getting 
ohmic dissipation, to establish the associated JFD relation and to explore the validity 
limits of the QCC principle. 

One may wonder what is the practical gain in achieving the above mentioned 
goals. Is it just a matter of doing 'mathematics' properly? A similar type of question 
is frequently asked with regard to the efforts to re-derive well known RMT results 
using semiclassical methods. The answer to such questions should be clear: It is 
not possible to analyze non-generic (or non-universal) features unless one possess a 
thorough understanding of a generic theory along with its limitations. In the future. 
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our intention is to analyze circumstances that go beyond these hmitations, and to look 
for genuine QM effects p3[ . 

We are using the term 'generic' frequently, and it is now appropriate to define what 
do we mean by that. The answer is as follows: In order to understand a phenomena 
(energy spreading and dissipation in the present case) it is a common practice to make 
the maximum simplifications possible. Then we get a theory with a minimal number 
of parameters. It turns out that the theory of energy spreading involves the minimal 
number of two dimensionless parameters. We are going to consider any additional 
(dimensionless) parameter as non-generic. Finally, it should be clear that some of 
our predictions concerning energy spreading are completely non-trivial. The kernel 
Pt(n\m), as well as its parametric version P{n\m), are accessible to numerical studies 
[ISL Em as well as to real experiments. 
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2. Energy surfaces and eigenstates 

We consider a system that is described by an Hamiltonian 7i((5, P; x) where (Q, P) are 
canonical variables and a; is a parameter. The phase space volume that corresponds 
to a the energy surface TL{Q^ P;x) — E is 

n = Q{E;x) = J ^^ OiE ~ n{Q, P; x)) (1) 

where d is the number of degrees of freedom. Measuring phase-space volume in units of 
{2itTiY is insignificant classically, but very convenient upon quantization. The density 
of phase space cells will be denoted by g{E) — d^Q.{E]x). The energy surface that 
corresponds to a phase space volume n will be denoted by \'n{x)) and its energy will 
be denoted by En{x). We assume a simple phase space topology such that for given 
n and x corresponds a unique energy surface. Thus 

\n{x)) ^ {{Q,P)\ n{Q,P;x)^E,,{x)} (2) 

The microcanonical distribution which is supported by \n{x)) is 

PnAQ^P) = -^MniQ,P;x)~E,,ix)) = S{n{HiQ,P;x))-n) (3) 

In the QM case the energy becomes quantized, and the mean level density is related 
to the classical density of phase space cells. By Weyl law we have: 

i ^ Y,S{E-E,,) = g{E) (4) 

n 

Thus, upon quantization, the variable n becomes a level-index, and Pn,x{Q, P) should 
be interpreted as the Wigner function that corresponds to the eigenstate \n{x)). With 
these definitions we will be able to address the QM theory and the classical theory 
simultaneously. We shall use from now on an admixture of classical and quantum- 
mechanical jargon. This should not cause any confusion. The QM discussion however 
is postponed to later sections. The following discussion is purely classical. 

Let us consider a set oi parametrically related energy surfaces \nix)) that enclose 
the same phase space volume n. By differentiation of the expression fl{E{x);x) = n 
with respect to the parameter x one obtains: 

SE = -Fix)Sx, F{x) EE /-^\ (5) 

The angular brackets denote microcanonical average over all the phase space points 
that satisfy 7i(Q, P; x) — E. Later we shall see that the quantity F{x) has the meaning 
of a generalized (conservative) force. Having F{x) — for any x is equivalent to having 
^{E; x) which is independent of x. Such is the case for a gas particle which is affected 
by collisions with a small rigid body that is being translated inside a large cavity. 
We shall refer to the latter example as the 'piston' example. See Fig.|l] and Sec.0 for 
more details. In order to simplify notations we shall assume, with almost no loss of 
generality, that indeed n{E;x) is independent of x. It is also useful to define 

-<*--' ^(-f)-(-f). 

We can define a parametric correlation scale Sx"^ that is associated with the function 
J-"((3, P; x). For the 'piston' example (to be discussed in Sec.^) it is just the penetration 
distance into the 'piston' upon collision (the effective 'thickness' of the wall). If either 
(Q, P) or x become time-dependent, then J^ becomes a fluctuating quantity. The 
nature of these fluctuations is discussed in the next section. 
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3. Fluctuations 

For a given x = x{t) and initial conditions ((5(0), P{0)) we may find the time-history 
{Q{t)P{t)). We shall also use the notations 

£{t) = n{Q{t),P{t);x{t)) 

T{t) = T{Q{t),P{t);x{t)) 

The correlator of the fluctuating force is 

C{t,T) ^ {Ht)Ht + r)) (7) 

The angular brackets denote microcanonical average over the initial (t=0) phase-space 
point ((3(0), P{0)). In the next two paragraphs we are going to discuss the statistical 
properties of of the fluctuating force J-'(i). First we consider the time independent 
case (x = const), and then the time dependent case. 

If a; = const, then £{t) = E is a, constant of the motion, and C{t,T) = Ce{t) is 
independent of t. It is assumed that the dynamics is chaotic, and consequently the 
the stochastic force looks like white noise. The fluctuations spectrum (7e(w) is defined 
as the Fourier transform of (7e(t). The intensity of the fluctuations is characterized 
by the parameter 

/•oo 

C^{T)dT = Ce(w=0) (8) 

The fluctuations are also characterized by a short correlation time t^, = 0^(0) / Ce{Q) ■ 
For generic Hamiltonian system it is natural to identify t„, with the ergodic time t^^^. 
However, in speciflc applications we may have r^i < t^^g- For the 'piston' example, 
that will be discussed in Sec.R, the correlation time t^j is equal to the duration of a 
collision with the wall, while i^rg is determined by the ballistic time r,,!. 

If X is time dependent rather than a constant, for example x{t) ~ Vt, then for any 
finite V and < i the actual distribution of {Q{t),P{t)) is no longer microcanonical. 
The statistical properties of the fluctuating force !F{t) are expected to be different 
from the V=0 case. The average {^{t)) is no longer expected to be zero. Rather, we 
shall argue (See (p2) that {J^it)) — —fiV, where fj, is the dissipation coefficient. This 
implies that the correlator C{t,T) acquires an offset {nV)^. The offset term can be 
neglected for a limited time t < t^^^ 

ifrc = '^/(m^)^ [classical breaktime], (9) 

provided {i^VY <C (7e(0). The latter condition, which is equivalent to having 

r^i ^ tf,c [non-trivial slowness condition]. (10) 

implies that the velocity V should be small enough. There is another possible reason 
for the correlator C{t, t) to be different from C^{t). Loss of correlation may be either 
due to the dynamics of {Q{t), P{t)) or else due to the parametric change of x{t). The 
correlation time which is associated with the dynamics is t^,. The correlation time 
which is associated with the parametric time dependence is r^' = Sx^^/V . We always 
assume that 

r^i ^ rf [trivial slowness condition] . (11) 

meaning that loss of correlations is predominantly determined by the chaotic nature of 
the dynamics rather than by the (slow) parametric change of the Hamiltonian. Since 
we always assume that that the classical slowness conditions ((lOh and dO)) are being 
satisfled, it follows that we can make the approximation 

C{t,T) = C{t) « Ce(t) fort<t,„ (12) 
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Figure 2. Schematic illustration of the dynamics. An initially localized 
distribution is launched in phase space. For a limited time (left plot) it travels 
upon the initial energy surface. But then it departs from it. After much longer 
times (right plot) the evolving distribution is concentrated across an instantaneous 
energy surface. 



4. Energy spreading and dissipation 



= -{F{x{t)) + T{t)) X 



For x(t) — const the energy £{t) is a constant of the motion. For time dependent x(t) 
and any particular time-iiistory we can write 

d£ _ m 

H ^ ~dt 

The first term implies reversible change of energy due to a conservative force that 
equals F(x). In what follows we shall see that the second term is responsible for an 
irreversible dissipation process. Integrating over time and using dq) we get 

£{t) = E{x{t)) - V I T{t)dt (13) 



where E{x) = Em{x) is the energy that correspond to the initial phase-space volume 
m. The difference E{x{t))—E{x{Q)) is due to the reversible work done by the 
generalized force F{x). If we disregard the fluctuating term, then we come to the 
conclusion that the trajectory is approximately bounded to the evolving energy surface 
\n{x{t))). Thus the phase-space volume n = n{£{t),x{t)) is an approximate constant 
of the motion, the so-called 'adiabatic invariant'. Using (|3[) we can estimate the 
energy dispersion which is associated with the fluctuating force: 



mt) ~ Eixit))^) = v dt' 

Jo 



t' 



C{t',T)dT 



(14) 



Hence we get a crossover from ballistic to diffusive behavior: 



mt) ~ E{x{t))f) . 


. C^{Q)-{Vtf 


for t < T,, 


(15) 


{{£{t) - E{x{t))f) . 


s 2D^t 


for T,i < t < tf„ 


(16) 



The ballistic spreading on short time scales just reflects the parametric change of the 
energy surfaces. This is the essence o f the sudden approximation which is illustrated 
in FigJ^ and further explained in App. Appendix A . The diffusive spreading on longer 
times reflects the deviation from the adiabatic approximation. See illustration in Fig.^ 
and further details in App. Appendix A[ The diffusion coefficient is 

rt 1 



^B - IV I C^{r)dT 



V' 



(17) 
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If Cjs{t) is short range in nature, then Z?e wiU tend eventuaUy to the weU defined 
constant value which is indicated in the right hand side of (p^). Note that the adiabatic 
approximation becomes exact in the formal limit V ^ 0, keeping Vt constant. 

For intermediate times {Q{t), P{t)) are distributed ergodically across the evolving 
energy surface, within a shell of thickness \/2D^t. We shall argue later (Sec.||) that 
more generally, for any t 3> torg, the spreading profile p{E) obeys the following diffusion 
equation 

For simplicity we assume here that there is no conservative work {F{x)=Q). The 
energy dependence of the diffusion process implies a systematic growth of the mean 
energy (£(t)) = J Ep{E)dE. Namely, 

Substituting D^ = I^'e^^ and integrating by parts one obtains Q — fJ,V^, along with 
J-I? relation that can be written schematically as /i = J^[i>]. The result for p depends 
on p{E). If p{E) is well concentrated around some energy E^ one obtains 

115 

Me = - „, -7r-={giE)v-E) fmicrocanonical version] (20) 

2g[E) dE 

Another, more familiar variation of the J-I? relation is obtained if one assumes a 
canonical distribution p{E) oc g{E)exp{~E/{kB,T)). By substitution into (p^QJ), or 
simply by canonical averaging over (£y), one obtains: 

/iT = — ; Vt [canonical version] (21) 

where i^^ is related to Ct^{t)^ and the latter is defined the same way as Ce(t), but 
with canonical rather than microcanonical averaging. 

Having energy dissipation implies that for V j^ the fluctuating quantity J-{t) 
has a non-zero average: From (|l3|), and recalling that we assume F{x) = 0, we get 
Q = {T{t))V. Therefore Q = ^V^ implies 

{Ht)) = -PV (22) 

In case of the 'piston' example V is the velocity of the 'piston' and {^{t)) is 
the associated friction force. In case of conductivity calculation (See Sec. 1.9), the 
parameter x represents (time-dependent) magnetic flux via a ring, V is the electro- 
motive-force, and {^{t)) is the drift velocity. 

5. Quantal energy-spreading: Linear response theory 

At first sight it seems that the classical derivation in SecJ4| applies also to the 
QM case provided !F{t) is treated as an operator. This is essentially the so-called 
'linear response theory'. The only approximation involved is C(t,T) « Ce(t). This 
approximation should be valid as long as the evolving state pt{Q,P) is similar to 
the initial microcanonical preparation. It is more difficult to satisfy this condition 
in the QM case. The similarity pt{E) w Po{E) is not a sufficient condition: It is 
also required that off-diagonal elements of the probability-matrix could be ignored, 
meaning that a superposition could be treated as if it were an incoherent mixture 
of the corresponding encrgy-eigenstates. The classical considerations (Sec|3|) lead to 
the time restriction t ^ ijrc- The QM considerations will lead to a stronger time 
restriction t <^ tp^^. Accordingly, the classical slowness condition (|l(]| ) is replaced by: 

Tai <C iprt [quantal slowness condition]. (23) 
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The determination of ip^t, which is the breaktime for QM perturbation theory, will 
be discussed in later sections. For non-slow velocities, (meaning throughout this 
paper that (M) and (O) are satisfied but (03) is violated), the following elementary 
considerations does not apply. The quantal slowness condition (p3) should not be 
confused with the QM-adiabaticity condition (to be discussed later). QM-adiabaticity 
requires extremely slow velocities. 

The QM version of the derivation in Sec.0 gives a classical look-alike result for the 
energy spreading. The only implicit modification is that the classical Ce(t) should be 
replaced by the corresponding QM object. For the purpose of concise presentation, 
the formula for the energy spreading can be written as follows: 

ft pt /•+00 I 

SE'' ^ V^ C^it2~ti)dtidt2 = Vh / —C^iio) Ftiiu) (24) 

Jo Jo J-oo ^TT 

where 

Ftiiv) = t-{smc{ujt/2)f (25) 

Now Ce(t) is a QM object, and its Fourier transform C^{lo) can be expressed as 



C^{^) 



I 



dn 

dx 



2nS[Lu- ^" /"' ) (26) 



One observes that the power-spectrum of the QM fluctuations has a discrete 
nature, and consequently the correlation function Ce(t) is characterized by the the 
additional time scale in- We assume that h is reasonably small such that r„, <^ in- 
Correspondence considerations imply that the quantal Ce(t) is similar to the classical 
Cb(t) as long as i ^ tn- Equivalently, as long as i <C tn the discrete nature of the 
quantal C'e(w) can be ignored, and we can effectively use the classical C'e(w). Recall 
that the power-spectrum of the classical fluctuations looks like that of white noise: 
It satisfies C'e(w) « i^e for \uj\ <C l/r^i and decays rapidly to zero outside of this 
regime. Thus, for t <C r^i we can make the replacement Ft{uj) — > t, and we obtain 
the ballistic result 6E^ = CE{0)-{Vt)^, while for t 3> r^ we can make the replacement 
Ft{uj) -^ 2tt6{uj), and we get then the diffusive behavior SE'^ — VeV^I. 

We can get a semiclassical estimate for the matrix elements in (Pq ) by exploiting 
the correspondence that has been mentioned above |Q. The function C'E(a;) is 
assumed to be vanishingly small for cj ^ 1/t„i which implies that {d'H/dx)nm is a 
banded matrix. Energy levels are coupled by matrix elements provided \En-~Em\ < Af, 
where 

Af, = -^ = band width (27) 

For LU <C l/r^i the smoothed C'E(a;) should be equal to the classical noise intensity fs- 
Consequently one obtains the following estimate for individual matrix elements within 
the band: 



a^ 



dn 

dx 



' A 

^ TTT"^ foi" \En-E„,\<Ah (28) 



It is important to specify the minimal number of (generic) parameters that are involved 
in the above analysis. In the classical problem there are just two generic parameters: 
Namely, r^i and i^E- Quantum-mechanics requires the specification of two additional 
parameters: Namely, the band width Ab and the mean level spacing A. The associated 
dimensionless parameters are: 

WpR = scaled velocity = y 2Z?e Td / At, (29) 

b = scaled band width = At / A (30) 
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The specification of A is not dynamically significant as long as t ^ t^. Longer 
times are required in order to resolve individual energy levels. Thus we come to the 
conclusion that in the time regime i <C tn there is a single generic dimensionless 
parameter, namely Vp^, that controls QCC. We shall see that the QM definition of 



slowness (23) can be cast into the form Upn <^ 1. 



6. Quantal energy spreading: The conventional FGR picture 

Equation ( p4|) for the QM energy spreading can be derived using linear response theory, 
i.e. by following the same steps as in Sec.^ However, the simplicity of linear response 
theory is lost once we try to formulate a controlled version of it. It is difficult to derive 
and to get a good understanding for the breaktime scale tp^. It is better to use the 
conventional version of timc-dcpcndcnt first-order perturbation theory (FOPT), and 
to view the energy spreading as arising from transitions between energy levels. 

The choice of basis for the representation of the dynamics is a crucial step in 
the analysis. The proper basis for the understanding of energy spreading is the x- 
dependent set of eigenstates \n{x)) of the Hamiltonian 7i((5, P; x{t)). This is the basis 
that we are going to use later in this paper. In a sense we are going to introduce an 
improved version of FGR picture. However, for sake of completeness, we would like 
to discuss in this section the capabilities and the limitations of the conventional FGR 
picture. The conventional FGR picture is using a fixed basis that is determined by the 
unperturbed Hamiltonian T-C{Q,P;x{0)). It should be clear that transitions between 
unperturbed energy levels reflect reduced-energy-chang es rather tha n actual-energy- 



changes (see corresponding classical definitions in App. Appendix A ). Therefore the 
description of the crossover from ballistic to diffusive energy spreading is out-of-reach 
for this version of perturbation theory. 

The conventional FGR picture can be used in order to determine the diffusion 
coefficient -De, as well as the perturbative breaktime T^^t. A detailed derivation 



can be found in App . Appendix E ) . We use the notation rprt rather than tprt, 
because a different version of perturbation theory is involved here. The final result is 
De = ^v^^V'^, with the effective noise intensity: 

Ce{t) F{t) S{t) dT (31) 

where F{t) is the correlation function of the driving source, S{t) = exp(— (r/2)t) is 
the survival amplitude, and Tp^ = 1/E. The introduction of S{t) is a common ad-hoc 



improvement of equation (B.3). Such type of improvement is used in other contexts 
to get the Wigner-Weisskopf Lorentzian line shape. It approximates the effect of 
higher orders of time-dependent perturbation theory. It is also important to specify 
the validity condition of the FGR picture. It is not difficult to be convinced that the 



requirement (B.4) can be relaxed slightly, and the actual condition is: 

FGR-condition: Either t„, or Tc <C Tp.t (32) 

Here Tq characterizes the correlation function F(t), namely, it is the correlation time 
of the driving source. It should be realized that having Tc ^ r^i constitutes an obvious 



variation of the trivial slowness condition (|ll|). Furthermore, using (B.5) one can easily 
conclude that a necessary condition for the applicability of FGR picture is Vp„ ^ 1, 
which coincides with the quantal slowness condition (p3). Thus, it is not possible in 
principle to apply the FGR picture in the limit ?i — > 0. 

The consistency of the FGR result (|3^) with the linear response result (|4|) is not 
obvious. It is true that ( p4|) gives a crossover from ballistic to diffusive behavior, where 
indeed D^ = \v'^*V'^ . But if one takes (|24| ) seriously for t ^ tn, one will come to the 
conclusion that this diffusion will stop due to recurrences [|2). The FGR result (pl|). 
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due to the presence of S{t), does not imply such a conclusion, provided t^^^ ^ in- A 
vanishingly small result for D^ is obtained only in the QM-adiabatic regime where we 
may have Tp^t ^ ^h- In the QM-adiabatic regime Landau-Zener transitions between 
neighboring levels become the predominant mechanism for energy spreading H, and 
the FGR picture becomes of minor importance. 

7. The 'piston' example and the v^rall formula 

The above picture and considerations become much more transparent once applied 
to cavities with moving walls. The Hamiltonian H — E{p) + V(x) describes the free 
motion of a 'gas particle', whose canonical coordinates are Q — x and P — p, inside a 
d-dimensional space which is confined by some boundary. Unless otherwise specified 
-^(p) = p^/2?7i where m is the mass of the gas particle. The corresponding velocity 
will be denoted by v^- The boundary is composed of wall-elements, and may have few 
components. For example it may consist of some 'static' component that defines an 
interior space, and an additional 'moving' component that defines an excluded space 
of an impenetrable 'piston' as in Fig.|l|. The displacement of the moving wall-elements 
is parameterized by x. The gas particle undergoes elastic collisions with the boundary. 
The ballistic time will be denoted by Tb,. It is determined by the collision rate with 



the walls. The derivation in App. Appendix F gives the result: 



;^^(E^('-'-)h5^<l""'l>''= ''=« 



The d-dependent geometrical factor (|cos6'|) is defined in App. Appendix D . For the 
purpose of 'ballistic-time' definition we should take the total Area of all the wall 
elements. For the purpose of calculating an effective collision rate the effective Area 
should be defined as in (p5|). For the 'piston' example the total effective Area of the 
moving-faces of the 'piston' may be much smaller compared with the total Area of the 
walls, and consequently the effective time between collisions will be much larger than 
the ballistic time. 

For the later QM considerations it is essential to consider 'soft walls'. For 
concreteness we may assume that the wall is realized by a constant force field /. 
If z is a coordinate perpendicular to a wall elements, than the potential barrier is 
V{z) — for 2 < and V(z) = f-z for z > up to some maximal value V^^u well 
inside the barrier. We assume that the energy E of the particle is much lower than 
V„aii, and therefore the latter energy scale should be of no significance. For strongly 
chaotic billiards successive collisions with the 'piston' are uncorrelated and therefore 
the classical correlation time r^i is equal to the collision time with the wall. Namely 

r,, = i2mv^)/f (34) 

Obviously, in the hard wall limit / ^ oo we have Td -^ 0. The displacement of the 
walls is parameterized by some parameter x. With each surface-element ds we can 
associate a normal unit vector n, and a 'propagation velocity' which will be denoted 
by V(s). The latter is simply the derivative of the wall-element displacement by the 
controlling parameter x. The effective moving-wall area is defined as follows: 

Area = (b{n-Y)'^ds (35) 

Due to ergodicity there are two different strategies that can be applied in order to 
calculate p^- One possibility is to average J^{t)J^{t+T) over (x, p), as implied by 
the definition (0), and then to integrate over t. Schematically the calculation goes 
as follows: (.F^) simply equal to /^ multiplied by the ratio between the collision- 
volume (Area X (weTci)) and the total Volume. Note that the latter ratio simply equals 
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Tci/T^ai- The noise intensity is obtained by further multiplication with r^i. The proper 



calculation should be done as in App. Appendix E. The other possibility to calculate 



Ue is to write J^(t) as a sum over short impulses and to perform the averaging in time 
domain. The magnitude of the impulses is 2mvE- The noise intensity is simply equal 
to the square of the impulses multiplied by the collision rate. The exact calculation 



is done in App. Appendix F. Both approaches give obviously the same result: 

Area 



1^^ = 2(1 cos en — mX (36) 

Volume 

This result is quite general, but it assumes that successive collisions with the 'piston' 

are uncorrelated. More generally, correlations between successive collisions should be 

taken into account. The derivation can be done by following an essentially identical 

computation by Koonin B and the result is cast into the form 

z/e - i" i"ds2dsi (n-V(s2)) i/(s2,si) (n-V(si)) (37) 

If one ignores correlations between successive collisions, then one obtains i^(s2, Si) oc 
(5(s2 — Si) and ( p7| ) reduces to (|3q). Using the JT> relation one obtains the following 
generalized 'wall formula': 

He = 2(1 cos e*!) — mvE (38) 

Volume 

The familiar d=3 version of the wall formula is obtained by substituting (| cos(0)|) = 
1/2. As in the the case of i/e we can try to derive this result using a simple-minded 



time-domain approach. See App. Appendix F. It turns out that only half of the 



correct result is obtained. Alternatively, the correct result (p8D can be obtained 



by extending the standard 'kinetic' derivation. See App. Appendix G, The kinetic 
derivation demonstrates that (pq) is more general than it seems at first sight. It 
applies to any velocity- momentum dispersion relation. The mass m = (dv/dp)^^ may 
be energy dependent. 

The QM calculation of z^°" requires the knowledge of the quantal Ce(t), and 
we should also have a proper understanding of F{t). For the time being let us 
assume that the effective Tc is much larger than r^,. It means that the transitions 
are resonance-limited, and detailed knowledge of F{t) becomes irrelevant. If the De- 
Broglie wavelength Ae = 2'Kh/{mvE) is much smaller compared with other (classical) 
scales, then it is expected to have QCC, as discussed previously with respect to (p6|). 



However, it is also possible to make a direct estimate of the matrix elements that 



appear in (|2q). See App. Appendix I, The results are in complete agreement with 



our semiclassical expectations. Indeed, the bandwidth is determined by the collision 



time with the (soft) walls (see App, Appendix J), and the power-law decay of matrix 



elements outside the band can be associated with the discontinuity (see App. Appendix 



E) in the derivative of the classical Ce(t) at r = 0. The expression for the effective 
noise intensity can be cast into the form of ( |37| ) with 

^(^^-^^) - "^^"^'^^3;P (Sincg|s2-s,|))'^ (39) 

The Sine function, as well as other notations, are defined in App. Appendix D| The 



QM result coincides with the classical result if Ae is small compared with the classical 
length scales that describe the s-dependence of n-V(s). It is also assumed that Ae 
is small compared with the surface radius-of-curvature, else further corrections are 
required [p|. 
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8. The route to stochastic behavior 

In this section we shall introduce a general phase-space formulation for the theory 
of energy-spreading. The main mathematical object of the study, namely the kernel 
Pt{n\m), will be defined. In order to go smoothly from the classical theory to the 
QM theory it is essential to use proper notations. From now on we use the variable 
n = i}{E) instead of E. See definitions in Sec.pl. The transition probability kernel 
Pt{n\m) is defined as the projection of an evolving state on the instantaneous set of 
energy-states. It is also possible to define a parametric kernel P{n\m). The latter 
depends on the displacement 6x but not on the actual time that it takes to realize 
this displacement. The definitions are: 

Pt{n\m) = trace( p^^xit) ^{t) Pni.,x(o) ) (40) 

P(n|m) = trace( Pn,x(t) Pm,x(o) ) (41) 

In the above definitions the initial energy-surface is \m{x{0))), and the associated 
phase-space density is Pm,x(o){Qj P)- In the QM-case |m(a;(0))) is an energy-eigenstate 
and Pm,a:(o)(Q7-F') is the associated Wigner function. The evolving surface/state is 
represented by U{t) Pm,x(o), where U{t) is either the classical Liouville propagator or 
its QM version. In the classical case it simply re-positions points in phase-space. In 
the QM case it propagates a Wigner function and it may have a more complicated 
structure. The trace operation is just a dQdP integral over phase-space. In the QM- 
case the definitions of P{n\m) and Pt{n\m) can be cast into a much simpler form using 
Dirac's notations: See ( p8| ) and (|59|). 

In the classical case the kernel P{n\m) reflects the parametric correlations between 
two sets of energy surfaces (Fig.^ upper left). Consequently non-Gaussian features 
may manifest themselves. An important special non-Gaussian feature is encountered 
in many specific examples where x affects only a tiny portion of the energy surface. 
(Fig.|| upper right). In the 'piston' example this is the case because {dTi/dx) — 
unless Q is near the face of the piston. Consequently P{n\m) will have a (^-singularity 
for n = m. 

The classical scenario for Pt{n\m) consists of three time regimes. For short times 
we have the classical sudden approximation: 

PtXn\m) « P{n\m) for t < r,i (42) 



See Fig.^ and App, Appendix A for more details. For longer times we have the classical 



adiabatic approximation, or more precisely we have diffusive spreading: 

Pt{n\m) Ri Gaussian(n— m) for t„,g ^ t -^ tf„ (43) 

For t 3> tfic the kernel Pt{n\m) is no-longer a narrow Gaussian that is centered around 
n = m. Using (||) it is easily observed that t > tt,^ is equivalent to Qt > \/2D^t, 
meaning that the systematic energy change becomes larger than the width of the 
spreading. Thus U^^ should be regarded as the breaktime for the classical adiabatic 
approximation. 

On time scales larger than ioig one may argue that the energy spreading is like a 
stochastic process, and consequently the diffusive behavior should persist beyond tfr^. 
A precise formulation of this point will be presented now. For t 3> tcrg the evolving 
surface U{t)\m{x{0))) becomes very convoluted due to 'mixing'. As long as one does 
not insist on looking for fine structures, pt{Q,P) can be replaced by its smeared 
version. Any 'tangential' non-homogeneity in phase space will be washed away due 
to the ergodic behavior, and therefore the smeared pt{Q,P) is fully characterized by 
the projected distribution pt{n). Obviously this statement is true only if i„g is much 
smaller than the time scale t{^^ that characterizes the 'transverse' spreading. Thus the 
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Figure 3. Upper Left: Phase space illustration of the initial and of the 
instantaneous set of parametric energy surfaces; Plot of the associated P(n|»n), 
where the classical behavior is indicate by the black lines, and the QM behavior is 
represented by the grey filling. Detailed QCC is assumed. In the QM case classical 
sharp-cutoffs are being smeared. Upper Right: Illustration of a typical non-generic 
feature. In the QM case the classical delta-singularity is being smeared. Lower 
Right: The same non-generic feature manifests itself in the 'piston' example. 
Lower Left: In the perturbative case there is no detailed QCC. The kernel is 
characterized by a core-tail structure. The tail is limited by the bandwidth of the 
coupling matrix-elements. If Sx is sufficiently small the core is just a kronecker's 
delta. 
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dynamics acquires the following stochastic property: 

Pt,+t2{n\m) K. j Pt^{n\n') dn' Pt,{n'\rn) provided ti,t2>i„g (44) 

Assuming that i„g ^ i&c it is possible to define an intermediate time ti = t/N, where 
N is some large integer, such that t„rg <C ii <C tfrc- Using the stochastic property 
(Q) one can write Pt{n\m) as a convolution of the N kernels Pt-^ (n|?7i). Applying the 
same considerations as in the derivation of the central limit theorem, we come to the 
conclusion that Pt{n\m) will become a spreading Gaussian that obeys the diffusion 
equation 

which is equivalent to (p^). Note that D^ — g{E)^D^. This description holds on 
time scales larger than i„g, irrespective of the detailed structure of Ptj(n|r?T,). Only 
the second moment of the latter is important for the determination of the diffusion 
coefficient. 

The argumentation in favor of long-time stochastic behavior is more subtle in the 
QM case. Using obvious notations the stochastic assumption ( p4| ) is: 

i|Ui,U,Jm)|') « ^(|(HU,,K)f)(|(^'|UtJ™)|') (46) 

n' 

Later we shall argue that U is a banded matrix. It is true in general that in the 
absence of correlations between successive unitary operations we will always have a 
stochastic diffusive behavior |23|. Thus, in order to establish a stochastic behavior in 
the QM case we should look for a time scale Tc that marks the loss of phase-correlation. 
For ii,i2 S> Tc we can argue that the off-diagonal 'interference' terms in the matrix 
multiplication Ui^Ufj will be averaged to zero. 

One way to establish the existence of a time scale t,. is just to assume irregular 
driving. As an example let us assume that we have a 'piston' that is pushed back and 
forth in arbitrary a-periodic fashion, meaning that x{t) becomes uncorrelated on a time 
scale that will be denoted by r^". The irregular driving is like noise and consequently 
the interference contribution is averaged to zero ||25| . For periodic driving, there may 
be limitation of diffusion due to 'localization' effect as in the quantum-kicked-rotator 
model [g^. The study of this latter issue is beyond the scope of this paper. As in 
the classical case we will be able to establish a diffusive behavior on an intermediate 
time scale. QM considerations will be limited either by a semiclassical breaktime t,^i 
or by a perturbative breaktime tp^t, which are analogous to the classical breaktime 
if,e. If we do not have the separation of time scales (r^'" ^ tp,t or t^" ^ t^^i) then we 
should wonder whether there is an intrinsic t^. An intrinsic t^ is expected to be either 
equal or larger than the classical time ie^g- Indeed, in the perturbative regime, where 
iorg ^ ^prt, it will be argued that effectively t^ ^ io,g- In the perturbative regime we 
are not able to give a general mathematical proof for having an effective Tc such that 
the separation of time scales requirement (tc <C iprt) is being satisfied. However, we 
are going to demonstrate that a crossover to a diffusive-growth of the second moment 
happens before the breaktime ip^ . The assumption that the diffusive behavior persists 
beyond ip^t with the same diffusion coefficient is the cornerstone of the common FGR 
picture. We are not going to study in this paper the general conditions for having 
stochastic-like behavior. 

The derivation of the classical J-I? relation consists of two steps: The first step 
establishes the local diffusive behavior for short (i <C itrc) time scales, and D^ 
is determined; The second step establishes the global stochastic behavior on large 
{t ^ ierg) time scales. The various time scales involved are illustrated in Fig.Q. The 
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Figure 4. Illustration of the various time scales involved in constructing either 
classical, semiclassical or perturbation theory of dissipation. We use the notation 
Tgc = ^^^sc/V. The accompanying table summarizes the associated requirements 
for the applicability of each of those theories. Note that the quantum mechanical 
definitions of slowness and of fastness are not complementary. Note also that 
slowness in the classical sense is always assumed. The rest of this paper (Sections 
9-20) is devoted to the study of the ballistic-diffusive crossover. This crossover is 
'captured' either by perturbation theory or by the semiclassical theory, provided 
the respective slowness or fastness conditions are satisfied. For simplicity we 
assume a generic system where Td can be identified with ierg ■ The persistence of 
the stochastic behavior for arbitrarily long times is assumed. We are not going 
to study in this paper the general conditions for having such long-time stochastic 
behavior. 



validity of the classical derivation depends on the slowness condition {W^ . The validity 
of the analogous QM theory is further restricted by the quantal slowness condition 
(P3|). However, an optional derivation of the JT) relation in the QM case can be based 
on semiclassical considerations. The limitations of the latter strategy are illustrated 
in Fig.y and further discussed in the next section. 



9. The semiclassical picture and detailed QCC 

The main objects of our discussion are the transition probability kernel Pt{n\m) and 
the parametric kernel P{n\m) which have been introduced in the previous section. 
Recall that we are measuring phase-space volume {n—Q{Ej} in units of {2nh)'^. This 
way we can obtain a 'classical approximation' for the QM kernel, simply by making n 
and m integer variables. If the 'classical approximation' is similar to the QM kernel, 
then we say that there is detailed QCC. If only the second-moment is similar, then we 
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say that there is restricted QCC. In the present section we are going to discuss the 
conditions for having detailed QCC, using simple semiclassical considerations. In the 
next paragraph we discuss the conditions for having detailed QCC in the computation 
of the parametric kernel P{n\m). Then we discuss the further restrictions on detailed 
QCC, that are associated with the computation of the actual kernel Pt{n\m). 

Wigner function pn,x(Q, P), unlike its classical microcanonical analog, has a non- 
trivial transverse structure. For a curved energy surface the transverse profile looks 
like Airy function and it is characterized by a width p3| 

1/3 

Asc - (£..(-) 1 (47) 




where e^i is a classical energy scale. For the 'piston' example e^i ^ E is the kinetic 
energy of the gas particle. The classical P{n\m) has a dispersion 



SE^, = 



\ 




(48) 



which characterizes the transverse distance between the intersecting energy-surfaces 
\m{x)) and \n{x+6x)). In the generic case, it should be legitimate to neglect the 
transverse profile of Wigner function provided SE^i 3> Age- This condition can be 
cast into the form 5x 3> Sx^c where 

Sxsc = -^^ oc h'^' (49) 

For the 'piston' example see [Q. Another important parametric scale is defined in a 
similar fashion: We shall see that it is not legitimate to ignore the transverse profile 
of Wigner function if SE^i < A(,. This latter condition can be cast into the form 
Sx <C (ScCprt where 

Jxprt = ^!L- = (50) 

Typically the two parametric scales are well separated (Sx^^t ^ Sxsc)- If we have 
dx <C feprt then the parametric kernel P(n\m) is characterized by a perturbative core- 
tail structure which is illustrated in Fig.|^ and further discussed in the next section. If 
we have Sx ^ Sxsc then the transverse profile of Wigner function can be ignored, and 
we get detailed QCC. Obviously, 'detailed QCC does not mean complete similarity. 
The classical kernel is typically characterized by various non-Gaussian features, such 
as sharp cutoffs, delta-singularities and cusps. These features are expected to be 
smeared in the QM case. The discussion of the latter issue is beyond the scope of the 
present paper jl^ . 

We turn now to discuss the actual transition probability kernel Pt{n\m). Here 
we encounter a new restriction on QCC: The evolving surface Li{t)\m) becomes more 
and more convoluted as a function of time. This is because of the mixing behavior 
that characterizes chaotic dynamics. For t ^ i^oi the intersections with a given 
instantaneous energy surface \n) become very dense, and associated QM features can 
no longer be ignored. The time scale i^ci can be related to the failure of the stationary 
phase approximation U%. 

The breaktime scale i^d of the semiclassical theory is analogous to the breaktime 
scale iprt of perturbation theory, as well as to the breaktime scale if,^ of the classical 
theory. In order to establish the crossover from ballistic to diffusive energy spreading 
using a semiclassical theory we should satisfy the condition r^ < t^^i. This velocity- 
independent condition is not very restrictive. On the other hand we should also 
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satisfy the condition 6x 3> <5xsci with Sx — Vt^i. The latter condition imphes that 
the appUcabihty of the semiclassical theory is restricted to relatively high velocities. 
We can define: 

Use = VDe T,i I Asc (51) 

If v^c ^ 1 then the above semiclassical analysis is applicable in order to analyze the 
crossover from ballistic to diffusive energy spreading. 

10. The perturbative picture and restricted QCC 

Detailed QCC between the quantal P(Ti|m) and the classical P(n\'m) is not guaranteed 
if bx < 5xsc- A-fortiori, this statement holds also for Pt{n\m). For sufficiently small 
parametric changes 5x, or for sufficiently short times t, perturbation theory becomes 
a useful tool for the analysis of these kernels. A detailed formulation of perturbation 
theory is postponed to later sections. In the present section we are going to sketch 
the main observations. We are going to argue that for small enough 5x there is no 
detailed QCC between the quantal and the classical kernels, but there is still restricted 
QCC that pertains to the second moment of the distributions. Large enough Sx 
is a necessary condition for getting detailed QCC. The following paragraph discuss 
the parametric evolution of P{n\m), and the rest of this section discuss the actual 
evolution of Pt{n\m). 

For extremely small Sx the parametric kernel P{n\m) has a standard 'first-order' 
perturbative structure, namely: 

P{n\m) « Sn,n + Tail(n-m) for Sx < Sx^"' (52) 

where (JxJ™ is defined as parametric change that is needed in order to mix neighboring 
levels. For larger values of Sx neighbor levels are mixed non-perturbatively and 
consequently we have a more complicated spreading profile: 

P{n\m) « Core(n— to) + Tail(7T,— m) for Sx <^ Sx^^^ (53) 

In the perturbative case {Sx ^ Sxp.t) the second moment of P{n\m) is generically 
determined by the 'tail'. It turns out that the QM expression for the second- moment 
is classical look-alike, and consequently restricted QCC is satisfied. The core of the 
quantal P(n\m) is of non-perturbative nature. The core is the component that is 
expected to become similar (eventually) to the classical P{n\m). A large perturbation 
Sx ^ (5xprt makes the core spill over the perturbative tail. If we have also Sx ^ Sxsc , 
then we can rely on detailed QCC in order to estimate P{n\m). The parametric scales 
Sx'^"' and Sxp^t are easily estimated in case of the 'piston' example. The displacement 
which is needed in order to mix levels is much smaller than De-Broglie wavelength, 
namely Jx^" ~ (A^+^/Area)^''^. The displacement which is needed in order to mix core 
and tail is much larger than De-Broglie wavelength, namely SXp^t ~ (tcoi/tci)^^^Ae. For 
a more careful discussion of these parametric scales see [Ol na| and the concluding 
section. 

The dynamical evolution of Pt{n\m) is related to the associated parametric 
evolution of P{n\m). We can define a perturbative time scale t^.^ which is analogous 
to feprt. For t ^ iprt the kernel Pt{n\m) is characterized by a core-tail structure 
that can be analyzed using perturbation theory. In particular we can determine the 
second moment of the energy distribution, and we can establish restricted QCC. If 
the second moment for the core-tail structure is proportional to t^, we shall say that 
there is a ballistic-like behavior. If it is proportional to t, we shall say that there is 
a diffusive-like behavior. In both cases the actual energy distribution is not classical- 
like, and therefore the term 'ballistic' and 'diffusive' should be used with care. We are 
going now to give a brief overview of the various scenarios in the time evolution of 
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Figure 5. The various crossovers in the time evolution of Pt{n\va). The vertical 
axis is x(t) = Vt. The parametric scales (Sxc™ and feprt are indicted by horizontal 
lines. The horizontal axis is the velocity Y . It is divided by vertical dashed lines 
to various velocity regimes. In each velocity regime there is a different dynamical 
route. The various crossovers are explained in the text and the various symbols 
are easily associated with having either Gaussian or some non-Gaussian spreading 
profile. In particular the perturbative spreading profile is either with or without 
non-trivial core, and its tail is either band-limited or resonance-limited. 



Pt(n|m). These are illustrated in Fig.||. In later sections we give a detailed account 
of the theory. 

For slow velocities such that r^i <C tp,t, there is a crossover from ballistic-like 
spreading to diffusive- like spreading at t ^ t^j. In spite of the lack of detailed QCC 
there is still restricted QCC as far as this ballistic-diffusive crossover is concerned. If 
the breakdown of perturbation theory happens before the Heisenberg time (iprt ^ ^h) 
it is implied that there is a second crossover at i ~ ip^t from a diffusive-like spreading 
to a genuine diffusive behavior. Once a stochastic behavior is established, the time 
scale tn for recurrences becomes non-effective, and we expect a long-time classical-like 
behavior. 

extremely slow velocities are defined by the the inequality tn ^ ipn- This 
inequality implies that there are QM recurrences before the expected crossover from 
diffusive-like spreading to genuine-diffusion. This is the QM adiabatic regime. In 
the t —^ oo limit Landau-Zener transitions will dominate the energy spreading, and 
consequently neither detailed nor restricted QCC is a-priori expected [g[. 

For fast velocities we have ip,t <?C r^, . There is a crossover at i ~ ip,t from ballistic- 
like spreading to a genuine ballistic behavior, and at i '^ t^j there is a second crossover 
from genuine-ballistic to genuine-diffusive spreading. The description of this classical- 
type crossover is out-of-reach for perturbation theory, but we can use the semiclassical 
picture instead. Note that the semiclassical definition of 'fastness' and the perturbative 
definition of 'slowness' imply that there is a 'gap' between the corresponding regimes. 
However, the interpolation is smooth, and therefore for simple systems surprises are 
not expected. 
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11. Actual versus Parametric Evolution 

The QM time evolution is governed by the time dependent Schroedinger equation 
with the time dependent Hamiltonian 7i(x(i)). In practice it is quite unnatural to 
use a fixed basis. For example, in case of the 'piston' example one may propose to 
use the fixed-basis that consists of the eigenf unctions of the empty cavity. However, 
the matrix elements of the 'piston' may be very large and even infinite if we assume 
impenetrable walls. Thus, it is much more natural to use the so called adiabatic 
basis, though the time evolution is not necessarily of adiabatic nature. The evolving 
state- vector is expanded as follows: 

m)) = ^a„(i)|n(x(i))) (54) 

n 

Using standard manipulation we obtain the Schroedinger-like equation 

—7- = "i'-^" "" ^ ¥ / y<lnrn{x{t)) a^, (55) 

m 

Were the off diagonal elements of W„m are 

(56) 



W„„i = -(n 



d \ hi 

■m 



dx 



dt I Er,-^ — Ej^ 

and we use the 'gauge' convention W„m=0 for n=m. (Only one parameter is being 
changed and therefore Berry's phase is not an issue). 

Equation (|55| ) will be now our starting point. It is defined in terms of W„m(a::) 
and in terms of a set of numbers {En}. Note that as long as t <C in we can ignore the 
dependence of En on the changing parameter x. The formal solution of ( p5|) will be 
written as follows: 

an{t) = ^U„„(t)a„(0) (57) 

m 

If all the En are set equal to the same constant, (or without loss of generality to 
zero), then ( pq ) describes the time evolution of a frozen wavefunction. In other words, 
(psh without the {£'„} is equivalent to the trivial equation dip/dt — 0. In this special 
case the formal solution ( |57| ) will be written with Tnm{x{t)) instead of U„„i(i). Note 
that if the {En} are taken away from (|55|), then x can be scaled out and therefore 
the X dependence rather than the t dependence become significant. The transition 
probability kernel and the parametric kernel can be written as: 

Ptin\m) = \Unm{t)\^ = \{n{x{t))\U{t)\m{x{0)))\^ (58) 

P{n\m) = |T„„(x)|2 = |(n(a:)|m(x(0)))p (59) 

From now on we shall refer to the i-dependent evolution which is represented by 
Unm(i) as the actual evolution (AE). To the i-dependent evolution which is represented 
by Tnm{x{t)) will shall refer as parametric evolution (PE). For PE the velocity x = V 
plays no role, and it can be scaled out from the above equation. Consequently, 
for PE, parametric scales and temporal scales are trivially related via the scaling 
transformation Sx — Vt. 

It is important to realize that in a certain sense, defined below, the AE coincided 
with the PE for short times t ^ t^^^. This is the QM-sudden approximation. The 
detailed picture is as follows: We start with some initial state |m). After time t 
there will be some non-vanishing probability to find the system in a certain energy 
range SE{t) around Em- As long as 6E{t) <C h/t the corresponding energy levels En 
within SE are not resolved. The latter condition defines a time interval t <C isdn- By 
definition, for t <C t^dn it is as if the energy-levels were degenerated. Therefore we can 
say that the AE coincides with the PE, implying that the evolving state (^^ remains 
approximately unchanged. The QM sudden approximation will be further discussed 
in Sec.^ 
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12. Application of perturbation theory 



We can use Equation (55) as a starting point for a standard first-order perturbation 
theory (FOPT). For short times, such that Pt{m\m) ^ 1, the transition probabihty 
from level m to level n is determined by the coupling strength |W„mP, by the energy 
difference {En— Em) and by the correlation function F{t). The latter describes loss 
of correlation between W„m(a:(0)) and W„m(a;(t)). It is defined via 

( W*„(t+r) W„,„(i) ) = \Wnmf Fir) (60) 

with the convention F{Q) = 1. Using FOPT one obtains the following result: 

2 



Pt(n|m) 



^ 

/ ** rt'. 



U. e'^V^* dt' 



/ / dt2dti F{t2-ti) e' 
Jo Jo 



(t2-tl) 

\ ^) Jo 
^ tPt ( ^"~^^"' ) X (^) for n ^ m (61) 

The function Ft{u!) describes the spectral content of the perturbation. For a constant 
perturbation {F{t) = 1) it is just given by equation (|25|). For a noisy perturbation 
F{tJ is characterized by some finite correlation-time Tc, and therefore the definition 
of Ft{uj) is modified as follows: 

^ ^ \ F{uj) for t > Tc ^ '^ 

where i^(w) is the Fourier transform of the correlation function F{t). Now it is a 
simple matter to calculate the second-moment of the spreading: 

/ + 00 7 
^C^cu) Ft{oj) (63) 

This result coincides with the linear-response result (|^ only if the coupling matrix- 
elements could have been treated as constant in time, meaning F(t) = 1 and accordingly 
Tc = oo. For t > Tc it becomes formally equivalent to the FGR result (|3l|). However, a 
practical equivalence seems unlikely because F{t) of (^0|) is not necessarily determined 
by the correlations of the external driving source. The critical discussion of this point 
is going to be the main issue of the subsequent sections. 

In order to use ( p3|) we should determine how F{t) look like, and in particular 
we should determine what is the correlation-time Tc. We postpone this discussion, 
and assume that F{t) and hence Tc are known from some calculation. The total 
transition probability is p{t) = ^^ P{n\m), where the prime indicates omission of the 
term n = m. FOPT is valid as long as p{t) <C 1. This defines a breaktime t'^^^ for 
the standard FOPT treatment. The above derivation imply that we can trust ( p3| ) 
only during the short time t <^ t\^. However, later we shall argue that with a proper 
(modified) definition of F{t) we can trust ( p3i ) during a longer time t <C iprt- The 
breaktime ip^t will be determined by using an improved perturbation theory (IMPT) . 

It is now possible to formulate the conditions for having restricted QCC. By 
'restricted' QCC we mean that only the second-moment of the spreading ( p3[ ) is being 
considered. It is essential to distinguish between two different possible scenarios: 

Resonance-limited transitions: Tc ^ r^i (64) 

Band-limited transitions: Tc <C t^ (65) 

For resonance-limited transitions, finite Tc has no consequence as far as SE^ is 
concerned. The crossover to diffusive behavior SE"^ ex t will happen at t ^ t„,. This 
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diffusive behavior will persist for t > Tc with the same diffusion coefficient. On the 
other hand, for band-limited transitions we will have at t ~ r^ a pre-mature crossover 
from ballistic to diffusive behavior. Consequently the classical result will be suppressed 
by a factor (tc/Tci) <C 1. This is due to the fact that the transitions between levels are 
limited not by the resonance width (embodied by F{u})), but rather by the band-width 
of the coupling matrix elements (embodied by Ce{lo)). 

We have realized that the perturbative result ( |63| ) can be used in order to establish 
a diffusive growth of the second moment. Obviously the applicability of this picture 
requires a separation of time scales: 

FGR-condition: Either r,,, or Tc <C tp^t (66) 

The long time stochastic behavior of the spreading is determined by the short-time 
dynamics, as explained in Section ^. The FGR condition guarantees that the diffusive 
growth of the second-moment is established before the breakdown of the short-time 
analysis. Therefore, the correct determination of the breaktime ip,t is extremely 
important, and it is going to be the main issue of the subsequent sections. 

13. The applicability regime of the standard FOPT treatment 

In order to have practical estimates for applicability regime of the standard FOPT 
treatment, we should look on the matrix W^^. This matrix is banded, and its 
elements satisfy: 



V Y 1 



SxT^ J {n^mY 



for |n-m| < b/2 (67) 



where fej'" — A/cr. From the above expression, once used in (61) for the calculation 
of the kernel P(n|rn), it follows that 5x'^ is the parametric change which is required 
in order to mix neighboring levels. Similarly, in the calculation of Pt{n\m) the related 
^qm _ 5xc/V is the time which is required in order to mix neighboring levels. Given 
two distant levels n and m, and taking the mixing on "small" scale into account, one 
realizes that Jo;)!'" also determines the correlation time Tc = r^" of the matrix-element 
'^nm{x{t)), as defined in (pQ). These observations can be summarized as follows: 

t'^^^ = Tc = T^"' for the standard FOPT. (68) 

The standard perturbative structure (^ of either P{n\m) or Pt{n\m) is maintained 
as long as neighboring levels are not being mixed. This structure obviously does 
not correspond to the classical structure since it is characterized by the non-classical 
energy scale A(,. Still, there is restricted QCC which is implied by (pSh. 

A sufficient condition for the applicability of the standard FOPT treatment is 
^RMT ^ 1- The argument goes as follows: By definition Wrmt ^ 1 implies t^, ^ t^"". 
Using (|6^) we observe that it is equivalent to r^i ^ t'^^^ . By definition iprt is either 



equal or larger than t'^^^. Therefore the FGR condition (66) is satisfied. The converse 
however is not true. Having f rmt S> 1 does not imply that the FGR condition cannot 
be satisfied. Therefore, for iJrmt S> 1, we cannot tell on the basis of the standard 
FOPT treatment whether or when there is a crossover to a diffusive behavior. We 
shall try to overcome this difficulty in the next sections. 

14. The over-simplified RMT (ORMT) picture 

Recall that the matrix W„m is a banded. The bandwidth 6=1,3... corresponds to 
diagonal, tridiagonal matrix and so on. In the spirit of RMT we can think of W„„i 
as a particular realization which is taken out from some large ensemble of (banded) 
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random matrices. In order to go beyond standard FOPT it is essential to further 
speciiy the cross-correlation between matrix elements. Following |10[ the simplest 
statistical assumption is absence of cross-correlations, namely, 

(W*,„,(i+r) W„™(t)) = if {n',m'}^{n,m} (69) 

Equation (|5g) with the statistical assumptions ( |69|) and (|60|), where Tc = r^", is a well 
defined RMT model. We shall refer to it as the over-simplified RMT (ORMT) picture. 
The main observation of ||l^ can be summarized as follows: The FOPT result (|63|), 
assuming an ORMT picture, can be trusted for classically long times. Namely, 

ip,, = tf„ and n = r^"" for the ORMT picture. (70) 

The ORMT picture reduces to FOPT and implies restricted QCC provided Wrmt <C 1. 
However, this condition is not satisfied in the classical limit {h ^ 0). In the classical 
limit Urmt S> 1, and consequently transitions are band-limited. In particular it follows 
from (pi^), that the classical diffusion (D^) is suppressed by a factor Tc/t^u leading to 

^r"" ~ ^^B for t,«MT » 1 (71) 

^RMT 

This result, which is the main result of |10|, is obviously inconceivable, because it is 
implied that the classical limit does not coincide with the classical result, and that 
the correspondence principle is actually violated. 

The ORMT picture predicts (for Wrmt ^ 1) a premature crossover from ballistic 
to diffusive behavior once Sx becomes larger than Sx"^"^. It is important to realize that 
([7l|), if it were true, would reflects a property of PE. This statement becomes more 
transparent if we write: 

= 2L>°™^ xt w ( ( ^ ) ) -^a^r X Sx (72) 




Exactly the same result would be obtained if we start with (p5|) without the the first 
term in the right hand side. The value of V has no significance in the above analysis. 

15. The core-tail structure 

There is no detailed QCC between the classical and the quantal Pt{n\m) for short 
times. We are going to explain that for a limited time Pt{n\m) consists of a core 
whose width will be denoted by 6(i), and a tail whose main component is contained 
within the bandwidth b. In order to analyze this core-tail structure we are going to 
use an improved perturbation theory (IMPT). The IMPT treatment assumes that out- 
of-band transitions {b/2 < |n— 77i|) can be neglected. The IMPT is useful as long as 
the second-moment of the energy distribution is predominated by the tail component. 
This determines the breaktime iprt. The breakdown of IMPT at t ~ ip^t happens 
once the core spills over the tail region, and the FOPT-like structure of Pt{n\m) is 
completely washed away. 

It is important to have a proper intuitive understanding of how the core is 
being formed. For t <C t^"" we have Pt{m\m) ~ 1 and Pt(n\m) <^ 1 for n ^ m. 
It means that the core width is b(t) — 1. For t ~ t^'" few levels are expected to 
be mixed by the perturbation, meaning that b{t) > 1. We may have the tendency 
to associate this mixing with an avoided crossing. However, this aspect should not 
be over-emphasized. The mixing of neighboring levels is not conditioned by having 
exceptionally small energy difference (having (Em+i—Em) <C A is not required). 
Moreover, one should not over-emphasize the importance of near-neighbor transitions, 
unless t^Lz 'C 1. (See further discussion of the QM-adiabatic regime in SecJ2y). If near- 
neighbor transitions were the dominant mechanism for energy spreading, it would 
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Figure 6. Schematic illustration of a generic core-tail spreading profile. The 
core-width b{t) is defined by the participation-ratio. The second-moment should 
satisfy b(t) <g; s{t) ^ b, where b is that bandwidth. In case of Pt{n\m) the 
tail becomes (for t > t^i) resonance-limited rather than band-limited. In the 
resonance-limited case the bandwidth b in the above figure should be replaced by 
{h/t)/A, and accordingly the requirement is b{t) <C s{t) <^ {h/t)/A. 



imply that b{t) « (t/r^''"). Rather, we shaU argue that the core develops much more 
rapidly, namely b{t) « (t/r^"")^- 

It is also important to have a proper intuitive understanding of how the tail is 
being formed. Let us assume for simplicity that only three levels to = 1, tti = 2 
and n — 100 are actually coupled by matrix elements. Let us start at i = with all 
the probability concentrated in ?n = 1. As we go at f = ii via an avoided-crossing 
of TO = 1 and TO = 2, the matrix element W„_i may change sign. This change of 
sign may be viewed as a loss of correlations. However, at the same time, assuming a 
diabatic transition, almost all the probability is transfered to to = 2 in such a way 
that W„2a2 at t > ii is strongly correlated with W„ la^ sd t < ti. It is implied 
that the effective correlation time Tj, for m -^ n transitions is larger than the time 
between avoided crossings. This is not captured by the statistical assumption ( |69| ) 
of the previous section. A proper transformation, that effectively removes the in-core 
transitions between TO-states, can be used for the purpose of tail-formation analysis. 
The associated perturbative treatment is characterized by an effective Tc ^ t^"" as 
well as by t'^^^ > t^"". 

For t > r^"" neighboring levels are being mixed and consequently the transition 
kernel acquires a non-trivial core-tail structure which is illustrated in Fig.pl The 
expression for Pt(n\m) can be written schematically as follows: 



Pt{n\m) « Core(n— ?7i) + Tail(7i— to) 
The kernel is characterized now by two scales: 



for t < i„ 



h{t) = core width = {^{Pt{n\m)Y \ 
s{t) = spreading = N [n—rnf' Pt{n\m) 



-1/2 



(73) 

(74) 
(75) 



such that 6(t) <C s{t) <C h. For t < r^" we have a trivial core with b{t) w 1, whereas 
for t 3> T^"' we have a non-trivial core with bit) ^ 1. The matrix elements satisfy 
(|W„mP) « l/{n—m)^. We shall see that in the 'band-limited tail' case we have 
Pt{n\m) ^ const/(n— to)^ up to the cutoff b, while for the 'resonance-limited tail' case 
we have Pt{n\m) ~ const/(n— to)^ up to the cutoff {fi/t)/ A. One should realize that 
the power-law behavior of the tail is 'fast' enough in order to guarantee that b{t) is 
independent of the tail's cutoff. The cutoff does not have any effect on the evolving 
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core. On the other hand, the second moment s(t), unhke b{t), is predominantly 
determined by the tail's cutoff, and it is independent of the core structure. 

16. An improved perturbation theory (IMPT) 

In order to extend perturbation theory beyond r^" it is essential to eliminate the non- 
perturbative transitions within the core. This can be done by making a transformation 
to an appropriate basis as follows: 

an{t) = ^T„™ c™(t) (76) 

m 

Tnrn = Tnm if \n-m\ < b' /2, else zero. (77) 

The amplitudes c„(t) satisfy the same Schroedinger equation as the a„(i), with a 
transformed matrix W. The general expression for W is 

W = t^ Wt - intt {dT/dt) + t^ET (78) 

where E is a diagonal matrix of the energies. This is a quite complicated expression. 
However, we are interested only in the core-to-tail transitions for which 

(W)„„ = (f twt)„™ for|n-m|>6' (79) 

(no approximation is involved). Once this transformation is performed the 'new' 
Schroedinger equation will be characterized by a new correlation time Tc and by a 
new perturbative time t\^^. Both Tc and t\^ will depend on the free parameter b' . Our 
choice of the course-graining parameter b' is not completely arbitrary. The restrictions 
are: 

• Unitarity is approximately preserved: b{t) <^ b' . 

• Core-to- Tail transitions are preserved: b' ^ b. 

• Long effective correlation time is attained: r^i <C Tc 

The feasibility of the last requirement deserves further discussion. Only matrix 
elements with |n— r7i| 3> b' are of interest, and therefore in the multiplication 
(Tt)„,„(W)„™(T)„™/ we can substitute (H) with {En-Em) replaced by {En'-E.^')- 
As b' becomes closer to 6, the matrix elements (W)n'm' become correlated on a time 
scale of the order r^' . This is because b' ~ b implies transformation to an x- independent 
basis. The time scale r^' has been defined in Sec||. As we change b' from b back to 
smaller values, we expect r^ to become smaller. By continuity, we expect no difficulty 
in satisftying the conditions b' <^ b and Tc 3> t^,, simultaneously. The validity of the 
improved perturbative treatment is further discussed in the next section. 

The usefuUness of the above transformation stems from the fact that due to the 
elimination of non-perturbative transitions within the core, t'^_^^ becomes much longer 
than T^"\ At the same time the information which is required in order to determine the 
second moment s{t) is not lost. We have |W„,„| « |W„,„| for core-to-tail transitions, 
and a practical approximation for the 'renormalized' spreading profile is 

«<"|™' - '-' + '« (^^) ^ {^)^ wftW '"" 

Breakdown of the improved perturbative treatment happens once the total transition 
probability becomes non-negligible (of order 1). Thus 

C = ib'r' X ^r (81) 

The behavior for |n— m| < b' is an artifact of the transformation and contains false 
information. However, for the calculation of the second moment only the tail is 
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significant. The tail is not affected by our transformation and therefore we will 
obtain the same result ( |63| ) for SE^ with one important modification: a different 
effective value for r^. Moreover, since b' is chosen such that t^ ^ r^i, it follows 
that the transitions are resonant-limited and consequently QCC is established also in 
the domain Urmt > 1- Obviously, at the same time we should satisfy the condition 
Tci <C iprt- It is easily verified that the latter condition cannot be satisfied if Upr > 1. 
This is not just a technical limitation of the IMPT strategy, but reflects a real difference 
between two distinct routes towards QCC. This point is further illuminated in the next 
section. 

17. Consequences of the IMPT treatment 

The IMPT is capable of giving information about the tail, and hence about the second 
moment. Given t, one wonders how much b' can be 'pushed down' without violating 
the validity conditions of our procedure. It is quite clear that b' ^ b{t) is a necessary 
condition for not having a breakdown of perturbation theory. If we assume that the 
energy-spreading-profile is characterized just by the single parameter b{t), then the 
condition b' ^ b{t) should be equivalent to t -^ t' ^^. Hence the following estimate is 
suggested: 

bit) - {t/rrr (82) 

We turn now to determine the Sx^^^ of the parametric evolution (PE) , and then 
the iprt of the actual evolution ( AE) . Recall that PE is obtained formally by ignoring 
the differences (En— Em), which implies that we can make in (pfl) the replacement 
Ft I— > t. Thus the tail of P{n\m) is band-limited and consequently the second moment 
is 

s(t)^ = 6 X (l/r^^")^ t^ [band-limited tail] (83) 

in agreement with the classical ballistic result (Iq). Our procedure for analyzing the 
core-tail structure of P{n\m) is meaningful as long as we have b{t) <C s{t) <C b. This 
defines an upper time limitation t^^^ which is related via Sx — Vt to the following 
parametric scale: 

Sx,,, = 6^/= Sxr - -^L= (84) 

At i = tprt we have b{t) ~ s{t) ~ 6, and we expect a crossover from a ballistic-like 
spreading to a genuine ballistic spreading. 

In the perturbative regime the AE departs from the PE once the energy scale Af, 
is resolved. This happens when t ^ r^i. The perturbative approach is applicable for 
the analysis of the crossover at i ~ r,,, provided Vt^i <^ feprt- This is precisely the 
condition Vp^ <^ 1. For t ^ r^, the tail becomes resonance limited (|n— to| < {h/t)/A) 
rather than band limited (|n— m| < b) and we obtain: 

s(t)^ = (1/r^'")^ in t [resonance-limited tail] (85) 

in agreement with the classical diffusive result (iq). Our procedure for analyzing the 
core-tail structure of Pt{n\m) is meaningful as long as we have b{t) ^ s{t) <C b. This 
defines a modified upper time limitation 

iprt = {rrT' t'J' = i—^j [applies to WpH < 1] (86) 

At t = ip,t we have b(t) ^ s{t) <C b, and we expect a crossover from a diffusive-like 
spreading to a genuine diffusive spreading. 
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18. Validity of the IMPT picture 

It is important to have a clear understanding of the difference between the IMPT 
picture and the ORMT picture. In both cases we can argue that P{n\m) has, for short 
times, a core-tail structure. The fundamental difference is the assumption concerning 
the effective Tc for core-to-tail transitions. The ORMT picture assumes that effectively 
Tc — T^", and equivalently the tails grow like 6x. This is the reason for having the pre- 
mature crossover to diffusive growth ( [7l|) of the second moment. The IMPT picture 
assumes that the effective Tf, is scale-dependent, and that the tails grow predominantly 
like 6x^. In other words, the tails grow as if we are still in the regime of FOPT. 
Therefore, from practical point of view, all we have to do in order to establish the 
validity of the IMPT picture is to verify that indeed the tails grow in a ballistic-like 
fashion (ex (5x^), and not in diffusive-like fashion (ex Sx). A related observation is that 
also the core- width b{t) grows as (ex Sx^). 

The argumentation (Sec. |l6|) in favor of the IMPT picture is not mathematically 
rigorous. It is therefore important to study specific examples. The obvious example to 
begin with has been defined by Wigner forty years ago ||l^, g^, |T^ . Let 7i = E + xB 
where E is a diagonal matrix and B is a banded random matrix. The IMPT picture 
should apply to the analysis of the PE of this Wigner model. Indeed, it is well known 
that P(n\m) for Wigner's model is a Lorentzian, and we may view this Lorentzian as 
a special case of core-tail structure. The width of the Wigner's Lorentzian, in energy 
units, is r = 2tt{{x-<j)/A)'^ x A. Thus we have indeed ex 6x^ for both the core and 
the tails. If the ORMT picture were true we would expect to get ex Sx dependence. 
All the results of the previous sections are consistent with the established results of 
Wigner. Having established the IMPT picture for PE, and observing that going from 
PE to AE involves no additional assumptions, it follows that we can safely proceed 
with the analysis as in Secp^. 

The validity of the the IMPT also has been verified numerically for the PE of 
the 'piston' example [^5|, and for the PE of a 2D nonlinear oscillator ||l^. It has 
been verified that indeed the tails grow in a ballistic- like fashion (ex Sx'^), and not 
in diffusive- like fashion (ex dx). Obviously, the assumption of having a structure- 
less core does not universally apply and also having b{t) ex 6x^ is a quite fragile 
result. If we want to have a better idea about the core structure we should apply, 
in any special example, specific (non-perturbative) considerations. In case of the 
'piston' example we can use semiclassical considerations |14[] in order to argue that 
the core has a Lorentzian shape whose width is H/t^^,. This semiclassical Lorentzian 
has nothing to do with Wigner's Lorentzian. The semiclassical Lorentzian is a purely 
non-perturbative structure. This structure is exposed provided (h/T^^^)/A <^ b{t), 
leading to the condition Sx :s> X^- Else we have a structure- less core whose width is 
characterized by the single parameter b{t) of (|8^). 

19. The quantum mechanical sudden approximation 

It is now appropriate to discuss the QM sudden approximation. For the perturbative 
scenario (wpR <C 1) we have already mentioned that the AE departs from the PE 
at tgdn = Tci, which is the time to resolve the energy scale A;,. In case of the non- 
perturbative scenario (vpu. ^1) there will be an earlier breakdown of the QM sudden 
approximation. This is because we have t^, 3> iprt and consequently at i = r^i we 
already have s{t) ^ b. Therefore i^dn should be defined as the time to resolve the 
energy scale which is associated with s{t). It leads to 

/ ^2 N 1/4 

i.a„ = b''\rrr^,y'' = f— ^j forz;p,»l (87) 



Theory of Quantum Dissipation 34 



perturbative route 


Kr < 1): 


tsdn ^ '^cl ^ t-prt 




At t = T„ 


b{t) < s(t) < 6 - (?iA)/A 


At t^t,,. 


6(t) - s(t) - (VO/A < & 


Non-perturbative 


route (wpR » 1): 


ip.t < i=d„ < T,, 




At i = tprt 


&(i)-s(i)-6<(n/t)/A 


At i = t,d„ 


& < s{t) - (?iA)/A 


At i = T,, 


&~(n/i)/A<s(i) 



Table 4. Various time scales in the route to stochastic behavior. 

The various time scales are summarized in Table 4. The non-perturbative crossover 
from genuine-ballistic to genuine-diffusive behavior in not trivial. If v^c ^ 1 we 
can relay on semiclassical considerations in order to establish the existence of this 
crossover. More generally, for v^j^ ^ 1, we would like to have an appropriate effective 
RMT model. This effective RMT model should support genuine-ballistic motion with 
an elastic scattering time r^i. 

20. The quantum mechanical adiabatic approximation 

The previous analysis has emphasized the role of core-to-tail transitions in energy 
spreading. An implicit assumption was that these transitions are not suppressed by 
recurrences. This is not true in the QM adiabatic regime {vi^z '^ !)■ Following M 
it is argued that energy spreading in the latter regime is dominated (eventually) by 
Landau-Zener transitions between near-neighbor levels. For completeness, the present 
section is devoted to the clarification of this observation. 

As a preliminary exercise it interesting to estimate the contribution of transitions 
between near-neighbor levels. The time scale that characterize these transitions is r^", 
and the 'step' size is A. Disregarding all other transitions, we have a random-walk 
process with diffusion coefficient (A)^/r^'", leading to 

Z?™ —D";^ [not appficable] (88) 

Thus for i>Lz S> 1 the contribution of near-neighbor (NN) transitions is indeed 
negligible as expected. In the QM adiabatic regime (ulz "C 1) the above result should 
be modified as follows 0: 

/ 1 \l-(/3/2) 

Dl^ « f— j i^^' for z;lz«1 (89) 

This result takes into account the no-trivial nature of Landau-Zener transitions and 
the statistics of the avoided-crossings. One should use /3 — 1 for the Gaussian unitary 
ensemble (GUE) and /3 = 2 for the Gaussian orthogonal ensemble (GOE). Recalling 
the stochastic considerations that lead to (p3) one deduces that the perturbative 
breaktime is 

/ 1 \/3/2 

tp.t = — rr (X T/-(i+(/3/2)) [applies to v^^^l] (90) 
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In the QM adiabatic regime energy spreading is dominated by near-neighbor level 
transitions for two distinct reasons. The first reason applies to the (3—1 case, namely 
D^^ » £)°'. The other reason is that D°' > D^'^''. In the latter inequality, L»^<^" 



is based on the FGR result (31). The FGR result becomes very small, compared 
with the classical result, once F{uj) becomes much narrower than the average level- 
spacing. The QM-adiabaticity condition t^ <C t^'" means that individual energy levels 
are being resolved before the breakdown of first-order perturbation theory. Having no 
'systematic' transitions to 'other' levels implies that the energy-distribution remains 
localized in the initial level for a very long time. The above argumentation implies 
that Z?^^ ^ D^"^^, meaning that for extremely slow velocities energy spreading, and 
the eventual breakdown of the QM adiabatic approximation, is predominantly due to 
Landau-Zener mechanism. 

21. Open questions and future directions 

The purpose of this paper was to make the first steps towards a theory for energy 
spreading and quantum dissipation. In particular we wanted to demonstrate that 
perturbation theory, and semiclassical theory have different regimes of validity. There 
are still a lot of open questions that have to be answered. 

An important issue is the specification of the general conditions for having a 
genuine stochastic behavior in the QM case. For fast velocities it is suggested (but 
not proved) that the stochastic behavior persists beyond the semiclassical breaktime. 
For slow velocities, it is suggested (but again not proved) that the stochastic behavior 
persists beyond the breaktime of perturbation theory. The latter suggestion is 
indirectly supported by common-wisdom and by various numerical experiments with 
banded matrices |^, |l^, Q. For the generic RMT picture, which is still lacking, it is 
implied that both, breakdown of perturbation theory and resolving the bandwidth of 
first-order transitions, are necessary conditions for having genuine stochastic behavior. 
In any case, stochasticity can be established if we assume irregular a-periodic driving 
with an appropriate correlation scale. For periodic driving, further considerations are 
required in order to analyze the possible manifestation of localization effect. 

A better understandingof the core-tail structure is required. Only in the case 
of Wigner's model ||l^, |l^, ^ we have an established result: Namely, the core-tail 
structure is simply a Lorentzian. For real systems the core-tail structure is not 
necessarily a Lorentzian |1J]. The determination of the border between the core 
and the tail may be problematic. One cannot exclude the existence of a distinct 
tail component, in the vicinity of the core, that does not grow like 6x^. A strongly 
related issue is to get an analytical understanding of the b' dependence of the effective 
('renormalized') correlation time Tc- 

The 'piston' example is non-generic in many respects. There are three classical 
length scales: The penetration distance upon collision with the piston; The mean 
path-length between collisions with the piston; And the ballistic length scale that 
characterizes the volume of the cavity. Quantum- mechanics adds two additional length 
scales: one is related to the Airy structure in the vicinity of the turning points, and 
the other is the De-Broglie wavelength. Having all these scales has some non-universal 
consequences |l4] that we have not considered in this paper. The application of the 
general theory of sections 8-20 to the 'piston' example is quite straightforward, but 
these non- universal features should be taken into account. In the generic theory there 
are only two parametric scales: The displacement dx"^"" that is needed in order to 
mix neighboring levels; And the displacement (JcCprt that is needed in order to mix the 
core with the tail. The former is much smaller than De-Broglie wavelength, and the 
latter is much larger than De-Broglie wavelength. It turns out that in the 'piston' 
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example there is a third, non-universal parametric scale iSx^u that roughly equals to 
De-Broglie wavelength Q. Consequently, the perturbative (slow- velocity) regime is 
further divided into a universal slow-velocity regime, and a non-universal slow-velocity 
regime. 

Of particular importance is the understanding of the hard-wall limit. In the 
generic theory t^ determines the bandwidth Ab of the matrix W„m. Having finite 
bandwidth is essential in order to understand that there is a crossover to a non- 
perturbative regime in the h —^ limit. We cannot treat W„„i as a banded matrix if we 
take first the limit t„i -^ 0. The walls of the 'piston' should be regarded as 'hard' once 
A(, becomes equal or larger than E. This is equivalent to having (classical) penetration 
distance smaller than De-Broglie wavelength. The consequence of taking the hard wall 
limit is that the non-perturbative regime {vp^ > 1) disappears. This state-of-affairs 
is possibly responsible to the illusion that a general theory for quantum dissipation 
can be base on a perturbative approach. At first sight it looks strange that hard-walls 
are 'better' for perturbation theory. It looks even more strange that for hard walls 
we cannot apply the semiclassical theory. In order to make the latter observation less 
strange recall that solving the one-dimensional Schroedinger equation near a sharp 
step, and then taking Ti ^ 0, never corresponds to the WKB approximation. 

The problem of quantum dissipation, in the sense of this paper, is a preliminary 
stage in the construction of a theory [Q of quantal Brownian motion (QBM). In the 
classical case it is known |1^ that the motion of a 'heavy' particle that is coupled 
to chaotic degrees-of-freedom is quite generally described by the classical Langevin 
equation. The effect of the environment is represented by a friction force plus a noise 
term. The friction leads to dissipation of energy and the noise term is essential for 
having diffusion. Furthermore, the friction coefficient is related to the noise intensity 
via the universal JT> relation. The fact that there is no general theory for quantum 
dissipation, and a-fortiori there is no general theory of QBM, has not been universally 
recognized in the literature. It is true that there is a vast literature that comes under 
those headings, but actually this literature is commonly based on an effective-bath 
approach. In previous studies [ p7| , Eq, E9| the common effective-bath strategy has been 
applied in order to develop a universal description of QBM and dephasing. Another 
possibility is to use an effective RMT bath pO] . The results of the latter study agree 
with |2^, |2^ . A future theory of QBM should clarify whether effective-bath methods 
universally apply. 
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Appendix A. The sudden and the adiabatic approximations 

For the time dependent Hamiltonian 'H{Q,P;x{t)) energy is not a constant of the 
motion. Changes in the actual energy £{t) reflect 'real' dynamical changes as well as 
parametric changes. Therefore it is useful to introduce the following definitions: 

£{t) = n{Q{t),p{ty,x{t)) (A.i) 

S'it) = niQ{t),P{t);xm (A.2) 

For simplicity we assume that the phase space volume 51(£'; x) is independent of x, 
thus F{x) = 0. The actual energy change can be calculated as follows: 

S£ = £{t) - £{0) ^ -V f T{t')dt' (A.3) 

Jo 

The actual energy change 5£ can be viewed as a sum of parametric-energy-change 5£o^ 

and reduced-energy-change 5£' . 

S£ = n{Q{t),P{t);x{t))-n{Q{0),P{0);x{0)) 
S£o = n{Q{t),P{t)-x{t))-n{Q{t),P{t)-xm 
S£' = H{Q{t),P{t);x{0))-n{Q{0),P{0);x{0)) 

The reduced energy change S£' reflects the deviation of (Q, P) from the original energy 
surface. It can be calculated as follows: 



5£' = F 



(A.4) 



T{t) xt - T{t')dt' 
Ja 

On the other hand, the actual energy change 5£ reflects the deviation of (Q, P) from 
the instantaneous energy surface. 

By inspection of the expressions for the reduced energy change we arrive at the 
conclusion that for short times we have the so called 'sudden approximation': 

5£' « fori<T,, (A.5) 

By inspection of the expression for the actual energy change we arrive at the conclusion 
that for long times we have the so called 'adiabatic approximation': 

^£ - fori<tf„ (A.6) 

The time evolution of an initially localized phase-space distribution pq(Q,P) is 
illustrated in Fig.H For short times we have in general non-stationary time evolution: 

Uit) poiQ, P) ^ poiQ, P) for t « r,i 

Here U{t) is the classical propagator of phase space points. However, if we operate 
with the same hi{t) on a microcanonical distribution, then 

U{t) Pe,x(o) {Q,P) ~ /Oe,x(o) (Q, P) for t < r,i 

Thus, the sudden approximation implies that for short times hi{t) can be replaced by 
unity if it operates on an initial microcanonical distribution. For long times we have 

Uit) PoiQ, P) - /5E,x(t)(g, P) for ie,, < t < it„ 

It is not required to start with a microcanonical distribution, unless r^i <C i ^ tcrg- 
The adiabatic approximation becomes worse and worse as time elapse due to the 
transverse spreading across the energy surface. t{„ is the breaktime of the adiabatic 
approximation. See discussion after (E^). 
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Appendix B. The conventional FGR Picture 

The simplest version of time-dependent perturbation theory is base on the 
approximated Hamiltonian 

n{Q,P]x{t)) « n{Q,P:x{0)) + ^x{t) (B.l) 

and using a fixed basis that is determined by the unperturbed Hamihonian 
T-C{Q,P;x{0)). For simphcity we set x(0)=0. A limitation that follows from using 
fixed basis is that the crossover from ballistic to diffusive spreading is out-of-reach for 
this version of perturbation theory. 

Another limitation of the present FGR picture stems from the fact that some 
additional assumptions should be imposed on x{t), else the treatment may be not 
valid. We shall assume that at any time x{t) <^ Sxl\ so that the expansion in 



(B.l) is valid. Moreover we shall assume that x is being changed in a an arbitrary 
a-periodic fashion, such that x(t) becomes uncorrelated on a time scale t^". The 
assumption r,,, ^ r^'" is implied by the trivial definition of slowness (|ll|) . The loss of 
velocity- velocity correlation will be described by a function F{t), with the convention 
F(0) = 1. This function should satisfy the normalization J F{T)dT — 0, else we 
will have an unbounded growth of {x{t) — a;(0))^. We still assume that the typical 
'velocity' is V, meaning that {i{t)x{t+T)) = V'^F{t). The above assumptions implies 
that the correlator {x{t)x{t+T)) = F(t) is well defined. Its Fourier transform satisfies 
F{llj) — V'^ F (uj) / ui'^ . Note that the normalization of the frequency distribution -F(w) 
is 1, while the normalization of F(ti;) is (Fr^")^. 

The FGR expression for the transition rate from the energy levels m to some 
other energy level n is: 

an' 

dx 

From here follows an expression for the diffusion constant. This expression is easily 
cast into a classical look-alike formula: 

11 1 f"" 

D^ = --^Y.'^nm{En-E„-,)^ ^ -V^ C^{T)F{T)dT (B.3) 



r - 1 

n 



'Er, Era , ^g 2) 



It is argued that (p3.2D is valid for any t > t^" provided there is a separation 
of time scales t^" ^ Tp^t. This is the 'golden- rule' condition, namely, breaktime 
of first-order perturbation theory should be after t^". We use the notation Tprt 
rather than tp^t in order to emphasize that fixed-basis perturbation theory is being 
used. The persistence of transitions with the same rate for t > Tp,t is guaranteed 
due the stochastic nature of the dynamics: The irregular driving is like noise, and 
interference contribution is averaged to zero once time intervals larger than t^'" are 
being composed. Let us assume further that Tp^t <C in, in order to guarantee that 
there are no recurrences irrespective of t^™. Gathering all our assumption together 
(classical slowness condition, FGR condition, and the non-recurrence condition) we 
get: 

T,i < r;?" < Tprt < in (B.4) 

It should be noted that the classical slowness condition implies resonance-limited 
transitions and hence restricted QCC is guaranteed. 

We turn now to get an actual estimate for the perturbative breaktime Tp^, and 
for the associated slowness condition. The total transition rate from a level m is 
T ^J2'n ^nm, and we have 

1 h'^ 

Tprt = f^ = — ,,n . ^ .y fixed-basis, -ypR < 1 (B.5) 
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The maximal value of Tp,t is attained if t^"' = r^j. The minimal value Tp^ — tprt 
is attained for t^" = tp^. One easily concludes that a necessary condition for the 
applicability of the FGR picture is 

V < / "' QM definition of slowness (B.6) 

The latter condition is always violated in the classical limit. Thus, it is not possible in 



principle to establish QCC in the limit ?i ^ by using FGR picture. Note that (B.6) 
is equivalent to (|2^) leading to Vp^ ^ 1- 



Appendix C. Alternative derivation of the 'sudden time' 

Simple considerations based on fixed-basis perturbation theory can be used in order to 
derive a result for t^^^. By 'sudden approximation' we mean that for a limited time we 
can 'ignore' the dynamical changes that are generated by the Hamiltonian. A necessary 
condition is that the transition probability between levels is much less than unity, or 
equivalently t <^ Tp^. It is essential to use the breaktime Tp^ oi fixed-basis perturbation 
theory, since we want to avoid the fake transitions due to non-trivial parametric- 
evolution. However, t <C Tp.t is not a sufficient condition. If |f/;) is a superposition 
of few levels, then we should also require that the corresponding energies will not be 
resolved. Namely we should have \En—Em\t <C h for any n and m in the superposition. 
We are interested in the way in which energy is re-distributed. Therefore, as long as 
first order perturbation theory applies, it is not important whether energy levels are 
resolved unless they are directly coupled {\En—Em\ < Ab). It follows that 

i,dn = minimum( rprt , t,, ) (C.l) 

where Tp^t is the breaktime of fixed-basis perturbation theory. It is a trivial matter 
to estimate the transition probabilities. Since we are interested in short times 



(i < T^i < T^") we can substitute x{t) — Vt in (B.l), and we get: 

for n ^ m (C.2) 



transition probability 



^Kf^ "'' 



Note that the FGR expression for the transition probability, namely Tnmt, is valid 
in a different time regime {t ^ t^" ^ t^i). The transition- probability to each of the 
levels within the band is approximately the same, and upon multiplication by &, an 
expression for the total transition-probability is easily obtained. It can be written as 
(t/Tprt)^ where 



fc2 \ 1/4 



v^y 



fixed-basis, Vpr 3> 1 (C3) 



If the velocity V is not slow (in the sense of (B_^)), then we have i,d„ = Tprt < Td. For 



slow velocities we have rp,t ^ t^ and therefore t^dn = Ta- In the latter case expression 



(C.3) underestimates the perturbative breaktime: Once t 3> t^i the transitions become 
resonance-limited, and consequently a diffusive-like spreading develops. Thus, in the 
slow velocity regime the FGR expression (B.5) for the perturbative breaktime Tp^t 



should be used. Taking into account the restrictions on r^" one observes that going 



from (C.2) to (B.J;) does not involve any discontinuity. 
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Appendix D. Spherical coordinates and related results 

The solid angle in d dimension and the volume element in spherical coordinates are 

^^ ^ rw) - 2, 2., 4., ... (D.l) 

dn^ r^'^dr = dna-i {sm{O)f-^d0 r'^-'^dr (D.2) 

The following results are easily derived: 

Note the relation 2(| cos(0)|> = (d+l)(| cos(6l)|3). 

It is useful to define a generalized cosine function by averaging exp(m-r) over the 
orientation of the unit vector n. The averaging is easily performed by using spherical 
coordinates: 

Cos(r) = (e^-=-(e)) = 2^-ir(f) %^ "" (D.5) 

It is also useful to define a generalized sine function as follows: 

Cos'fr) 1 d Jd(r) 

Sinc(r) = -^^2m = L 2lr(|+l) ^ (D.6) 

The function (Sinc)^ will appear in an integral over a (d—l) dimensional surface. It 
will be possible to replace it by an effective delta-function: 

(SincWs.-s,|))^ ^ {^f^'-^^^^^Sis^-s,) (D.7) 

In order to derive this result one should use the following: 



CO 



^ (Sinc(r))2 n,_, r'-'dr = ^ -^ (r(|))2 (D.9) 



^2-1 ^ ^2^ 
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Appendix E. Collisions with a wall: phase space approach 

The position of a particle inside a cavity can be described by Q = (z,x_l) where z 
is perpendicular to the surface. The potential which is experience by the particle is 
assumed to be V{Q) = inside the cavity where z < 0, and V{Q) = f-z for z > 0. 
Let the variable x parameterize the perpendicular displacement of the surface (hence 
n-V = 1). With this parameterization we we have simply !F{Q,P;x) — f provided 
z > 0, else T{Q, P; x) = 0. For an isolated collision we have 

z(t±ir) = z±^.(lr)-i^(irf (E.l) 

The correlation function is 

P Jdp^Jdp^JdK^JdzS(E^(^^+f-z)) 
C^ir) = / .\ — (E-2) 

The integration in the numerator is restricted by the conditions z(t + -^t) > and 
z{t — ir) > 0. The dx± integration gives an Area factor. After the dz integration one 
is left with a dp integration that represents the volume p^ + 2(i/r)p2 + (i/r)^ < 2mE. 
Using spherical coordinates one obtains 

, . _ Are^.■ f ■ 2 J fld-iisiney-^de J p''-^ dp ^ ^ ( ^ 
^^^' ' Volume- if^rf {{2mEy'^r / E " "^^^ [r^, 

The d6 and the dp integrations in the numerator are restricted by the conditions 
{I fr) / {2mEy^^ < coa0 < 1 and (^fr)/ cosO < p < {2TnEy'^ respectively. The 
noise intensity is 

ly^ ^ 2{\cos0\'') -pp^ m^vl (E.3) 

Volume 

The properly normalized Ce(t) equals zero for |r| > 1 and otherwise can be expressed 
using an hypergeometrical function as follows: 



C4r) 



^" JtKcosB 
d+\ 1 I \ ' 

2d \t^ ) 


19 

2f 


\ cos 9 J 

^ /d~l d d+ 
^ V 2 ' 2' 2 


for d = 2 
for d = 3 


-kr 


\ J (^arctan (^^J ^ - 
I |(|rp- 




3|, 


-1 + 2) 





(E.4) 



I t is much easier to obtain an explicit expression for Ct(t). The J(energy) function in 
( |E.2| ) should be replaced by exp(— energy /(fee T)). The integration is factorized, and 
one obtains 

The noise intensity v^^ can be found either by integration over this function, or else by 
performing thermal average over Vj^. In both cases the result is: 

^^ ^ 4 Are^ ^^^ {2mk^Ty/^ (E.6) 

^y^T Volume 
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Appendix F. Collisions with a wall: time domain approach 

Consider the motion of a particle inside a cliaotic cavity, we can define and estimate 
the collision rate with a wall element as follows: 

Above teal is the time of of a collision, and ds is the area of the wall element. The 
derivation of last equality is based on ergodicity considerations which we are going to 
explain now. Let us define the coordinate z as in the previous appendix. We have 

where t^j correspond to the crossing times of the z = surface, and t^^i — (^ii+^^O/S 
is the time of the collision. The duration of the collision Ai^„i = (i+j — t~j) is assumed 
to be extremely short. Equivalently we can write ^ 5(i— icoi) = |we | cos0| 5{z{t)). In 



order to get (F.l), the latter expression should be averaged over the (implicit) initial 
conditions ((3(0), P(0)). As always (...) indicates this type of averaging. Due to the 
chaos, the coordinate z and cos 6* can be treated as independent variables. Due to 
ergodicity we have {d{z{t))) = {6{z))e = ds/Volume. As always (...)e stands for a 



microcanonical average. For (|cos6'|) one may substitute (D.3). 

Exactly the same procedure can be used in order to estimate ly^- The first step is 
to realize that dH/dx = (n-V) x / for z > 0, and zero otherwise. The fluctuating force 
J^{t) = —dH/dx looks like a train of short impulses. The duration of each impulse is 
At^ai — '^'n^\{vz)coi\/ f- Therefore we can write 

^(i) = -^ = Y. [ ("■'^) 2™«- ^o^(^-) ] ^(^-^-0 (F-2) 

col 

The intensity of the fluctuations due to this random-like sequence of delta-impulses is 
j^E = (^ [ (n-V) 2mvE cos(0,„,) ]^ <5(t-t„,) ) 



The last expression is manipulated as in ( FA ) and leads to the correct result 



It is more illuminating to repeat the last derivation using a kinetic point of 
view. The velocity of the particle prior to the collision with the wall element is v„„i. 
The orientation of the wall element is represented by a normal unit vector n. The 
velocity of the wall element is denoted by V. Only the n-V component of the wall 
velocity is significant. After the collision the z component of the particle's velocity is 
n-(— v^„i -I- 2V). The corresponding energy gain due to the collision is 

A^,„, = -2m (n-V) (n-(v,„, - V)) . (F.3) 

The stochastic-like force T{t) is a sum over impulses that are created by collisions 
with the surface of the cavity. The energy gain due to this collisions is obtained by 
integrating —Vy^{t) over time (see Sec. 4). It follows that each collision involves an 
impulse —AE^ai/V leading to 

Ht) = Et 2"^ ("-V) (n-(vco> - V)) ] 6{t - t.) (F.4) 

col 

For 1^=0 we can assume that the velocities w^oi are uncorrelated. For a moving slab- 
shaped 'piston' half of the collisions will be from the 'left side' with n-V = —V, and 
half of the collisions will be from the 'right' side with n-V = +V. Consequently we 
have {T{t)) = 0. More generally, we have {T{t)) — for any parametric deformation 
that keeps constant the total volume of the cavity. 
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It is now tempting to make an attempt for a direct estimate of {^{t)) for F 7^ 0. 
Let us assume for a moment that we have a slab-shaped 'piston' which is moving from 
'left' to 'right'. At first sight it seems that friction is due to the fact that energy-gain 
due to collisions from the 'right' is larger than energy-loss due to collisions from the 



'left'. Our "hypothesis" implies that the 0(y) term in (FA) should be responsible for 



the friction. However, it turns out that only half of the correct result is obtained: 
/^2m(n-V)2<5(i-t,„,)\ = {\ cos e\) ^^- mv^ = ^^^^ 



We conclude that the first term in (F.4), the term that includes v^^i, has a non-zero 
average. Once 1^ 7^ we can no-longer treat the collisions as uncorrelated. The 
statistical effect (that has been emphasized in the derivation of the J^ relation) is 
important also in time-domain analysis. 

Thus we see that the !FD indirect approach is a quite powerful tool. This becomes 
manifestly evident once we consider a variation of the above example: If successive 
collisions with the 'piston' are correlated, for example due to bouncing behavior, then 
it is still a relatively easy task to estimate Vj^ for the V=0 case, and then to obtain 
He via the JT> relation. On the other hand, a direct evaluation of fi-E using kinetic 
considerations is extremely difficult , because in calculating {J-'{t)) it is essential to 
take into account subtle correlations between successive collisions. 

Appendix G. The wall formula: the standard kinetic approach 

Here we generalize the derivation of the wall formula to any dimension, and to 
any dispersion relation E = E{p). We use the standard kinetic approach. The 
momentum probability distribution of the particles inside the cavity is uniform in space 
and p(p) in momentum. (Later we are using many-particles jargon, but the actual 
meaning is always probabilistic). We assume isotropic distribution, and therefore 
the energy distribution is just p{E) — g{E)p{p). The velocity which is associated 
with a momentum p is v=dE/dp. The mass is defined in the differential sense 
m{E) = [dv/dp)-^. 

With a surface area ds we associate a normal unit vector n, and a displacement 
velocity V. The particles that actually collide with the wall-element during a time 
interval dt must satisfy the condition (v— V)-n > 0. Those with a velocity v are 
contained in a volume element dVolumc = ((v— V)-n) dtds. The energy gain due to a 



collision, assuming V <^ v, and expanding up to 0{V'^), is still given by (F.3). The 
total energy change is 

dE = p{p)dp X dVolume x A£',„i (G.l) 

The above integral can be manipulated as follows: We make a transformation to an 
integration variable p' = p— ttiV and correspondingly v' = v— V. Then we expand 
p{p'+m'V) = p(p') + [d p{p') / dE)nw' -Y , and throw away the first term since it is 
associated with reversible work. From here on we omit the primes. Due to symmetry 
consideration we may replace v-V by (v-n)(V-n). Now we obtain 

dE = - I ^^dpm^{^V■nf\^^■■a\^ dsdt (G.2) 

The latter integral is easily cast into the form of (O) with the identification 

De = ^i^m^.^ kv.nfds (G.3) 

Volume J 

The latter expression incorporates integration over all the wall-elements. 
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Appendix H. Chaotic eigenstates 

The average density of energy eigenstates for a particle in d dimensional cavity is 

1^1 }^ n, k^'^ (H.i) 

where A; is the wavenumber and v^ is the corresponding velocity of the particle. An 
eigenstate that corresponds to a chaotic cavity looks like a 'random wave'. More 
precisely, it has the same statistical properties as those of an un-constrained random 
superposition of plane waves whose wavenumber is |k| = /c. It is characterized by the 
correlation function 

(^k(xi)i^h(x2)) = — ^Cos(fc|x2-xi|) (H.2) 

Volume 



The function Cos is defined in App. Appendix D. We can constrain a random wave 



to be zero along a boundary by making an antisymmetric superposition of an un- 
constrained random wave with its mirror image. The normal derivative along the 
boundary ip{s) = n-Vf/" satisfies 

2fc2 

(^h(si)(Pr(s2)) = — Sinc(/i:|s2-si|) (H.3) 

Volume 

This expression ignores curvature effect. We assume that the boundary's radius-of- 

curvature is very large compared with De-Broglie wavelength, as well as compared 

with the distance of interest r = |s2 — Si|. 

Appendix I. Matrix elements for hard v^ralls 

The position of a particle in the vicinity of a wall element can be described by 
Q = (z, s) where s is a surface coordinate and z is a perpendicular 'radial' coordinate. 
The potential which is experience by the particle is assumed to be V{Q) = inside 
the cavity where z < 0, and V{Q) = f-z for z > up to some maximal value Vwaii 
well inside the barrier. We assume E <C V„aii and therefore V^^n should have no 
significance. The limit f ^ oo corresponds to hard walls. For deformation field V(s), 
defined such that V(5x is the displacement of a wall element, one obtains 

^ = -(n(s).V(s)) V„.„ Siz) (I.l) 

The orientation of each wall element is represented by the unit vector n(s). The 
logarithmic derivative of the wavefunction along the boundary is ip{s)/^(s) where 
(p{s) = n-Vtfj. For z > the wavefunction ipiQ) is a decaying exponential. Hence 
the logarithmic derivative of the wavefunction on the boundary should be equal to 
~y/2rnV^i/h. Consequently one obtains the following expression for the matrix 
elements that are associated with the deformation: 

— j = -— j) ipn{s)(pm{s) (n-'V)ds (1.2) 



For uncorrelated chaotic eigenstates we get by squaring (L2) and using ( [H.3| ), the 
following expression: 

2\ / 1 X 2 //t,_\2\ 2 



m 

dx 



f I \ f(hkf 



V Volume / \ m 



X ® ® Sinc(fc„|s2-si|) Sinc(A:,„|s2-Si|) (n-V(si))(isi (n-V(s2))ds2 

If De-Broglie wavelength is a small scale, then the difference \kn — km\ will have no 
significance in the calculation. Soft walls are essential in order to have finite bandwidth 
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(see next appendix). A classical look-alike result for i/^ — C'e(O) is obtained by setting 
k,i "^ kjn ^ k and then multiplying by 27r?i/A. Assuming that Ae is a small scale, it is 
possible to approximate the (Sinc)^ by a delta function (see (D.7)) and consequently 
the semiclassical estimate that follows from (Eq) for individual matrix elements is 
recovered: 



dn 

dx 



1 



{hkf /2(|cos(6')|s 



Volume 



^d 



Area K 



(1.3) 



Compared with non-chaotic eigenstate, for which Sine h^ 1, there is a factor l/vA/", 
where M = (Area/A^^^) is the number of correlated regions on the surface of the 
cavity. 



Appendix J. Finite bandwidth due to soft walls 

Equation (O) is valid for hard walls. For soft walls \{dT-L/dx)nm\'^ should be multiplied 
by a suppression factor Cgoft < 1- This suppression factor can be estimated by 
considering a plane wave incident upon the wall B . The result is 

(J.l) 



COS( 



COS( 



where e = {k/{2mf /Tl^Y'^Y is the scaled energy of the eigenstate. The reflection 
phase-shift is obtained via tan(0) = {A/A')^/e. We define A = Airy{—e), and A' and 
A" are the respective first and second derivatives. If </)„ is very close to (j)m then we 
can make the approximation 

2 

(J.2) 



{A')^ - {AA" 



[A') 

In the hard wall limit e ^ 0, 
unity as expected. If e„ and e^ 
uncorrelated and we obtain 

1 



'+e{AY 
and A" ■ 



0, consequently this factor will be equal 



are not very close, then </>„ and (pm can be treated as 



1 



1 



The bandwidth can be estimated by finding the difference £"„ — E,n for which the last 
expression becomes of order unity. This gives the result A;, = h/r^i where t^, is the 
collision time with the wall. This is in agreement with our semiclassical expectations. 
We learn that outside the band, the matrix element have a power-law decay. This 
result looks of semiclassical nature since it adds an h independent factor to Ce(w). A 
l/w^ tail of Ce(w) should be associated with a discontinuity of the first derivative of 



Ce(t) at r = 0. See App. Appendix E 
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